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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 


Non-SI  units  of  measurement  used 
(metric)  units  as  follows: 

In  this  report  can  be 

converted  to  SI 

Multlolv 

Bv 

To  Obtain 

feet 

0.3048 

meters 

Inches 

2.54 

centimeters 

miles 

1.609 

kilometers 
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SBEACH:  NUMERICAL  MODEL  FOR  STMUt ATINC  STORM- INDUCED  BEACH  CHANGE 
REPORT  2:  NUMERICAL  FORMULATION  AND  MODEL  TESTS 

PART  I :  INTRODUCTION 
Problem  Statement  and 

1.  Beach  stabilization  and  coastal  flood  protection  are  major  areas  of 
concern  In  the  field  of  coastal  engineering.  Beach  erosion,  accretion,  and 
change  In  the  offshore  bottom  topography  occur  naturally  through  the  transport 
of  sediment  by  waves  and  currents .  Additional  changes  result  from  perturba¬ 
tions  Introduced  by  coastal  structures,  beach  fills,  and  other  engineering 
activities.-1  Beach  change  Is  controlled  by  wind,  waves,  currents,  water  level, 
sediment  characteristics  and  supply,  and  constraints  on  sediment  movement  such 
as  those  imposed  by  coastal  structures.  These  controlling  factors  are 
nonlinear  and  have  great  variability  In  space  and  time.  Although  it  Is  a 
challenging  problem  to  predict  the  course  of  beach  change,  such  estimations 
are  necessary  to  design  and  maintain  shore  protection  projects. 

2.  The  beach  can  be  considered  as  a  flexible  structure  that  provides 
protection  to  life  and  resources  along  the  coast.  However,  over  an  interval 
of  Just  a  few  hours,  storms  can  produce  serious  damage  and  life-threatening 
situations  by  rapid  erosion  of  the  beach  and  upland  inundation.  Quantitative 
estimation  of  storm- Induced  beach  change  (e.g.,  contour  retreat  and  eroded 
volume)  as  a  function  of  the  time-dependent  forcing  conditions  Is  therefore  of 
great  Importance.  Another  important  aspect  of  beach  response  to  storm  action 
is  the  recovery  phase  normally  occurring  at  the  end  of  a  storm  and  during 
periods  of  mild  wave  conditions.  As  wave  height  decreases  during  the  waning 
stages  of  a  storm,  vaves  transform  from  erosive  to  accretive,  and  the  fore¬ 
shore  Is  rebuilt.  It  Is  vital  to  describe  the  process  of  recovery  if  storm 
Impact  Is  to  be  assessed  and  the  amount  of  material  lost  from  the  beach 
quantified.  The  recovery  process  must  also  be  modeled  If  the  Impacts  of 
successive  storms  are  to  be  estimated. 

3.  Seasonal  changes  In  profile  shape  are  readily  apparent,  produced  by 
smaller,  longer  period  waves  in  the  summer  and  higher,  shorter  period  vaves  in 
the  winter.  Waves  associated  with  summer  conditions  tend  to  produce  sand 
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movement  directed  toward  shore  Co  build  up  the  subaerial  beach,  whereas  winter 
waves  typically  erode  the  Inshore  part  of  the  beach  profile  and  deposit 
material  offshore,  creating  one  or  more  longshore  bars.  It  Is  Important  to 
reproduce  these  seasonal  changes  In  profile  shape  if  predictions  or  simula¬ 
tions  are  to  be  valid  for  long  time  periods.  Also,  the  capability  to  estimate 
variability  In  shoreline  position  due  to  seasonal  and  shorter-tlme  changes  In 
wave  climate  provides  a  basis  for  managing  coastal  resources  rationally  and 
Identifying  If  and  when  a  net  loss  In  beach  volume  has  taken  place. 

4.  Initial  adjustments  of  beach  fill  to  wave  action  and  the  response  of 
beach  fills  to  storm  processes  are  of  central  Interest  in  the  design  of  fill 
schemes  for  artificial  nourishment.  Prediction  of  protective  dune  and  beach 
profile  evolution  for  different  geometric  fill  configurations  and  grain  sizes 
can  provide  guidance  for  an  efficient  benefit-cost  analysis  of  proposed 
nourishment  designs.  Different  nourishment  templates  must  be  evaluated  and 
compared  objectively  to  assist  in  finding  the  optimal  design  to  satisfy 
criteria  based  on  storm  protection  and  recreation  benefits. 

5.  Structures  located  in  the  nearshore  zone  can  alter  beach  evolution, 
restricting  sand  transport  In  a  manner  dependent  on  the  type,  size,  and 
location  of  the  structure.  In  particular,  potential  profile  change  on  a 
seawall -fronting  beach  must  be  estimated  to  determine  If  the  wall  will  be 
damaged  during  storm  attack.  Likewise,  the  Influence  of  the  wall  on  modifying 
the  beach  Is  of  Interest  in  project  design.  Knowledge  of  beach  change  Induced 
by  seawalls  will  help  in  their  design  for  optimal  performance  of  the  structure 
and  to  appraise  Its  Impact  on  the  beach  (cf.  Kraus  1988). 

6.  Prediction  of  beach  evolution  with  numerical  models  has  proven  to  be 
a  powerful  technique  that  can  be  applied  to  assist  In  determining  project 
design.  Models  provide  a  framework  for  developing  project  problem  formulation 
and  solution  statements,  organizing  data  collection  and  analysis,  evaluating 
alternative  designs,  and  optimizing  the  selected  design.  Problems  related  to 
beach  profile  response  during  storms,  initial  adjustments  of  a  beach  fill, 
seasonal  changes  In  profile  shape,  and  the  influence  of  seawalls  on  profile 
evolution  are  conveniently  evaluated  using  numerical  models. 

7.  This  report  Is  the  second  In  a  series  describing  development  of  a 
numerical  simulation  model  for  engineering  use  for  predicting  beach  profile 
evolution.  Including  the  formation  and  movement  of  major  morphologic  features 
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such  as  longshore  bars,  troughs,  and  berms.  The  model  Is  called  SBEACH,  an 
acronym  which  stands  for  £torm- induced  BEAch  EHange •  The  first  report  In  the 
SBEACH  series  (Larson  and  Kraus  1989a),  hereafter  called  Report  1,  described 
model  development,  including  the  approach,  empirical  foundation  of  the  model, 
and  initial  tests  and  verification  with  large-scale  laboratory  data  and 
limited  field  data.  The  objectives  of  the  present  report  are  to  describe  the 
numerical  solution  procedure  and  to  demonstrate  applicability  of  the  model 
with  newly  acquired  field  data  pertaining  to  realistic  situations  encountered 
in  engineering  practice.  Report  3  in  this  series  will  provide  details  on  the 
operation  of  SBEACH  in  the  form  of  a  User's  Manual. 

Review  of  Methodology 

8.  The  model  SBEACH  was  developed  using  data  obtained  from  experiments 
performed  in  large  wave  tanks  with  monochromatic  waves  of  prototype -scale 
heights  and  periods.  Report  1  gives  a  thorough  discussion  of  the  data  base, 
including  a  wide-ranging  empirical  analysis  of  the  data  to  establish  cauae- 
and-effect  relationships  between  the  waves  and  beach  profile  change.  Based  on 
this  analysis  and  physical  considerations ,  a  numerical  model  was  developed  and 
tested  for  simulating  profile  evolution  generated  in  the  large  wave  tanka. 

The  model  was  then  tested  using  field  data  obtained  at  CERC’s  Field  Research 
Facility  (FRF)  in  Duck,  North  Carolina,  and  was  shown  to  give  reasonable 
agreement  between  calculated  and  measured  bar  formation  and  movement. 

9.  In  this  report,  the  capability  of  the  model  to  predict  berm  and  dune 
erosion  is  evaluated  using  recently  acquired  field  data  from  both  the  US  East 
and  West  coasts.  Hypothetical  cases  are  also  simulated  to  deisonstrate 
applicability  and  potential  uses  of  the  model  to  predict  the  initial  adjust¬ 
ment  of  a  beach  fill  and  its  response  to  storm  action,  including  post-storm 
recovery.  Two  different  synthetic  storms  are  developed,  corresponding  to 
tropical  (hurricane)  and  extratropicsl  (northeaster)  events. 

10.  Improvements  in  the  numerical  model  have  been  made  beyond  those 
presented  in  Report  1.  The  wave  transformation  model  has  been  generalized  by 
using  complete  linear  wave  theory  everywhere  along  the  profile  without  any 
shallow  water  wave  approximations.  This  capability  becomes  important  in 
applications  encountered  in  engineering  practice.  Also,  It  is  now  possible  to 
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specify  an  arbitrary  angle  between  the  Incident  waves  and  the  bottom  contours, 
thus  taking  into  account  wave  refraction.  An  option  has  been  Included  to 
randomize  the  wave  input  to  better  represent  forcing  conditions  In  the  field. 
This  is  done  by  randomly  changing  the  input  wave  height  by  a  pre 'determined 
percentage  around  Its  Input  value.  Such  a  modification  of  the  wave  height 
makes  the  break  point  move  back  and  forth  and  thus  produces  smoother  profile 
features  more  in  accordance  with  the  topography  observed  In  the  field.  Wave 
runup  is  presumably  better  represented  in  this  process,  thereby  taking  another 
step  toward  improved  representation  of  sediment  movement  on  the  foreshore. 

Scope  of  Report 

11.  Part  I  provides  an  introduction  and  short  review  of  Report  1. 

Part  II  reviews  the  empirical  foundation  of  the  model,  with  information  on  the 
model's  capabilities  and  limitations,  as  well  as  related  subjects  for  future 
work.  Part  III  describes  the  wave  model  which  provides  input  to  the  empiri¬ 
cally-developed  transport  formulas  and  a  realistic  estimate  of  the  water  level 
under  storm  surge,  tide,  and  waves,  as  described  in  Part  II.  Part  IV  describes 
the  profile  change  model,  focusing  on  the  numerical  solution  scheme  and  data 
input  requirements.  Part  V  first  reviews  sensitivity  tests  of  the  model  and 
then  describes  two  case  studies,  one  for  the  US  East  coast  and  one  for  the  US 
West  coast.  In  Part  VI,  the  model  is  applied  for  hypothetical  examples  to 
demonstrate  a  range  of  uses  in  beach  fill  design.  Part  VII  provides  concluding 
observations  of  the  present  capabilities  and  limitations  of  SBEACH,  leading  to 
recommendations  for  further  study  and  development  of  the  model. 
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PART  II:  REVIEW  OF  EMPIRICAL  FOUMDATION 


12.  This  chapter  gives  a  summary  of  the  empirical  foundation  of  the 
beach  profile  change  numerical  simulation  model  S BEACH .  The  data  base  used  In 
model  development  Is  outlined,  and  the  relevance  of  laboratory  data  Is  dis¬ 
cussed  with  reference  to  field  conditions.  Salient  aspects  of  cross-shore 
sand  transport  characteristics  derived  from  the  data  are  presented,  focusing 
on  direction  and  distribution  of  transport.  Transport  regions  used  In  the 
model  are  Introduced  and  related  to  findings  In  nearshore  wave  dynamics.  In 
addition,  the  semi-empirical  cross-shore  sand  transport  rate  equations  are 
summarized. 


Data  Base 

13.  The  complexity  and  randomness  of  nearshore  sediment  processes,  the 
limited  capability  of  quantifying  nearshore  fluid  and  sediment  movement,  and 
the  limited  availability  of  input  data  in  realistic  applications  all  indicate 
that  numerical  models  for  simulating  beach  profile  change  in  engineering 
situations  must  be  based  on  empirical  relations.  In  Report  1,  it  was  argued 
that  a  macroscale  approach  of  model  development  based  on  a  firm  empirical 
foundation  provided  a  reasonable  approach  which  paralleled  the  regularities  in 
morphologic  change  exhibited  by  beaches.  These  changes  are  on  a  spatial  scale 
of  meters  and  on  a  time  scale  of  hours  to  days .  Such  regularity  contrasts 
markedly  from  the  randomness  and  complexity  of  individual  particle  movement. 
However,  the  lack  of  detailed  field  data  encompassing  high-resolution  measure¬ 
ments  in  time  and  space  together  vith  associated  vave  and  water  level  condi¬ 
tions  made  it  necessary  to  use  laboratory  data  to  enhance  development  of  such 
a  model. 

14.  Use  of  laboratory  data  must  address  the  problem  of  scaling  in 
attempts  to  transfer  observations  and  results  to  the  prototype.  It  is  desir¬ 
able  to  perform  experiments  at  a  scale  vhlch  is  close  to  that  of  waves 
producing  significant  change  on  a  natural  beach,  without  distorting  the  kine¬ 
matic  and  dynamic  interaction  between  fluid  and  sediment  motion.  However, 
only  a  small  number  of  experiment  programs  have  been  carried  out  on  a  scale 
comparable  to  the  field  situation  (Savllle  1957,  Kajima  et  al.  1982,  Velllnga 
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1982,  Dette  and  Uliczka  1987).  Data  from  the  first  two  laboratory  programs 
were  used  for  numerical  model  development  In  this  study  (Report  1)  because 
they  provided  an  extensive  range  of  conditions  that  reasonably  replicated  a 
naturally  sloping  beach.  In  the  combined  data  set,  a  large  number  of  wave  and 
beach  conditions  were  available  for  analysis,  and  the  empirical  relations 
obtained  were  less  likely  to  be  influenced  by  experiment-specific  conditions. 

15.  Model  development  was  largely  based  on  data  obtained  in  the 
prototype -scale  experiments,  whereas  model  verification  also  used  data  from  a 
natural  beach.  Larson  (1988)  and  Larson  and  Kraus  (1989a)  give  a  detailed 
discussion  of  the  two  laboratory  experiment  programs  together  with  the  results 
of  extensive  data  analysis  of  the  geometric  characteristics  of  the  profile  and 
their  relation  to  wave  and  sand-beach  properties. 

16.  In  this  report,  only  a  short  review  of  the  experiments  is  presented 
to  provide  a  background  for  discussion  of  the  numerical  model.  The  first  data 
set  was  obtained  in  experiments  performed  by  the  US  Army  Corps  of  Engineers 
(denoted  by  the  abbreviation  CE)  in  the  years  1956-57  and  1962  (Seville  1957, 
Caldwell  1959,  Kraus  and  Larson  1988)  at  Dalecarlla  Reservation,  Washington, 
DC.  The  second  data  set  was  derived  from  experiments  carried  out  at  the 
Central  Research  Institute  of  Electric  Power  Industry  (denoted  by  the  abbrevi¬ 
ation  CRIEPI)  in  Japan  (Kajima  et  al.  1982,  1983). 

CE  experiments 

17.  The  CE  experiments  were  conducted  in  a  concrete  tank  193.5  ■  long, 
4.6  m  wide,  and  6.1  m  deep  and  having  a  standard  operating  depth  of  4.6  m  (for 
details  of  the  experimental  facility  see  Raynor  and  Simmons  (1964)  and  Kraus 
and  Larson  (1988)).  Most  cases  were  started  from  a  plane  initial  beach  slope 
of  1:15,  and  waves  of  constant  height  and  period  were  applied  until  near- 
equlllbrlum  conditions  were  attained.  A  run  typically  lasted  40  to  60  hr,  and 
the  water  level  was  held  constant  for  all  cases  except  one  in  which  a  step- 
type  sinusoidal  variation  was  employed.  On  the  order  of  10  to  15  profile 
surveys  were  made  during  a  typical  case,  and  the  survey  interval  was  1.2  m. 

18.  Two  different  sands  were  used  with  median  grain  sizes  of  0.22  and 
0.40  mm.  Eighteen  distinct  cases  have  been  documented  (Kraus  and  Larson  1988) 
of  which  16  started  from  a  plane  slope.  Wave  height  varied  between  0.55  and 
1.68  a  in  the  horizontal  portion  of  the  tank,  and  wave  period  ranged  between 
3.75  and  16.0  sec.  Some  observations  were  made  of  the  breaking  wave  height 
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and  location  during  the  runs,  and  a  limited  number  of  wave  runup  measurements 
were  performed  as  well. 

CRIEPI  experiments 

19.  The  CR1EPI  experiments  were  similar  to  the  CE  experiments  in 
design,  and  wave  tank  dimensions  (205  m  x  3.4  m  x  6  m)  were  comparable  to 
those  of  the  CE  tank.  The  experiment  program  consisted  of  24  cates  with  a 
varying  Initial  beach  profile,  including  both  plane  Initial  slopes  and 
profiles  of  arbitrary  shape.  The  generated  wave  height  ranged  between  0.3  and 
1.1  ■  and  the  wave  period  between  3  and  12  sec.  All  cases  were  run  with  mono¬ 
chromatic  waves  and  a  fixed  water  level.  For  cases  starting  from  a  plane 
beach  slope  (17  of  the  24  cases),  the  initial  slope  varied  from  1:50  to  1:10. 
Two  different  sands  with  median  diameters  of  0.27  and  0.47  mm  were  used. 

20.  Wave  height  and  wave  setup  along  the  profile  were  measured.  The 
number  of  profile  surveys  made  during  a  run  was  somewhat  fewer  than  for  the  CE 
experiments,  but  the  spatial  resolution  of  the  measurements  was  considerably 
higher.  Refer  to  Report  1  for  details  of  data  analysis. 

Relevance  of  large  wave  tank  experiments 

21.  Although  scale  effects  In  the  large  wave  tanks  were  probably  small, 
some  differences  In  response  of  the  beach  profile  as  compared  to  a  natural 
beach  are  expected.  The  main  questions  concern  the  relevance  of  the  labora¬ 
tory  data;  that  Is,  have  the  main  mechanisms  governing  cross-shore  sand 
transport  on  a  natural  beach  been  reproduced  In  the  wave  tank,  and  have 
additional  effects  been  Induced  In  the  tank  environment  which  are  not  found  on 
a  beach?  One  procedure  for  evaluating  the  appropriateness  of  wave  tank  data 
Is  to  compare  observed  profile  evolution  with  that  found  In  the  field;  such 
comparisons  are  summarized  below. 

22.  The  speed  of  bar  movement  calculated  from  the  large  wave  tank  data 
(Report  1)  was  very  similar  to  observations  made  by  Sallenger,  Holman,  and 
Blrkemeler  (1985)  on  a  natural  beach  during  a  storm  event.  In  general,  the 
time  scale  of  profile  change,  both  onshore  and  offshore,  was  rapid  for  the 
field  measurements  and  in  accordance  with  what  was  observed  in  the  large  wave 
tank  experiments.  In  particular,  onshore  bar  movement  occurring  during  the 
post-storm  recovery  stage  took  place  rapidly  In  the  field  Just  as  was  noted 
during  the  accretionary  cases  In  the  large  wave  tanks. 
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23.  The  ratio  between  depth  to  trough  bottom  and  depth  to  bar  crest 
obtained  In  the  large  tanks  displayed  properties  in  good  agreement  with  field 
observations  reported  by  Shepard  (1950).  Bars  in  the  field  exposed  to  waves 
with  longer  wave  periods  showed  a  smaller  value  of  the  ratio  between  depth  to 
trough  and  depth  to  crest  in  accordance  with  the  wave  tank  data.  The  value  of 
this  ratio  was  preserved  during  offshore  bar  movement  in  the  large  wave  tank 
experiments  which  Sallenger,  Holman,  and  Blrkemeler  (1985)  also  found  to  be 
the  case  during  a  storm  event  recorded  in  the  field. 

24.  A  criterion  to  predict  overall  profile  response  (on-  or  offshore 
sand  movement)  based  on  the  sand  fall  speed  and  the  deepwater  wave  steepness 
was  derived  from  the  large  wave  tank  data.  This  criterion  proved  to  be  valid 
for  natural  beaches  if  the  mean  wave  height  was  used  in  the  criterion  as  the 
representative  statistical  wave  height.  The  two  empirical  coefficients  in  the 
formula  expressing  the  criterion  were  the  same  for  the  wave  tank  data  and 
field  data,  further  emphasizing  the  relevance  of  wave  tank  experiments  as  well 
as  the  apparent  lack  of  scale  effects  in  the  experiments. 

25.  Hallermeier  (1987)  pointed  out  the  pertinence  of  results  recorded 

'  during  large-scale  laboratory  tests.  He  found  quantitative  agreement  In  terms 
of  shoreline  retreat  and  erosive  geometry  between  a  large  wave  tank  case  and 
measured  erosion  during  a  storm  at  Bethany  Beach,  Delaware,  where  the  condi¬ 
tions  in  the  tank  and  in  the  field  were  similar.  The  capability  of  S BEACH  to 
predict  eroded  volume  as  veil  as  breakpoint  bar  movement  in  field  altuatlona 
also  supports  the  validity  of  large  wave  tank  studies  (Report  1) . 

26.  An  obvious  restriction  of  the  wave  tank  data,  which  is  an  advantage 
from  the  viewpoint  of  clarifying  cross -shore  transport  characteristics,  is  the 
absence  of  longshore  fluid  and  sand  motion.  However,  one  of  the  main  reasons 
for  conducting  profile  change  experiments  in  wave  tanks  is  to  gain  insight 
into  the  physical  processes  of  cross-shore  sand  movement.  The  large  vave  tank 
experiments  were  conducted  using  waves  with  a  constant  height  and  constant 
period  (monochromatic  waves)  during  a  run.  Thus,  it  is  an  underlying  assump¬ 
tion  that  quantitative  understanding  of  profile  evolution  may  be  obtained  by 
studying  the  response  to  monochromatic  waves  of  uniform  height.  In  the  model, 
the  effect  of  Irregular  waves  is  simulated  by  using  a  time  series  of  waves  of 
different  heights  and  periods .  This  procedure  is  probably  adequate  for 
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simulating  profile  change  by  individual  irregular  waves  but  neglects  the 
contribution  of  other  types  of  waves  such  as  infragravity  waves. 

27.  An  alternative  approach  is  to  employ  a  statistical  or  probabilistic 
description  of  the  incident  waves,  which  would  certainly  lead  to,  for  example, 
more  gently  sloping  bars  that  develop  in  the  field  as  opposed  to  the  steeper 
bars  that  develop  under  monochromatic  waves  of  constant  height  in  tanks.  A 
proper  limit  of  such  models  is  replication  of  profile  change  and  bar  growth 
under  an  extremely  narrow-banded  wave  height  and  period  distribution.  In 
development  of  SBEACH  there  was  concern  that  an  a  priori  probabilistic  or 
statistical  description  would  Introduce  smoothing  and  mask  or  obliterate  the 
significance  of  individual  wave  components.  Therefore,  the  approach  used  Is 
founded  on  profile  change  produced  under  elemental  monochromatic  waves  of 
constant  height  and  proceeds  to  describe  change  under  Irregular  waves  by 
superposition. 

28.  Randomly  varying  wave  height,  period,  and  direction  produce  a  break 
point  which  continually  changes  position,  leading  to  smoother  profile  develop¬ 
ment  with  features  less  pronounced  than  In  the  monochromatic  wave  case 
(Keulegan  1948)  .  In  the  large  wave  tank  experiments  beach  slopes  tended  to 
become  steep  as  the  profile  approached  equilibrium,  especially  around  the  bar 
or  at  the  foreshore  step  (Report  1) .  This  was  attributed  to  the  applied 
constant  wave  height,  wave  period,  and  water  level  which  promoted  gradual 
break  point  movement  and  the  consequent  build-up  of  steep  beach  slopes.  In 
summary,  a  fundamental  assumption  guiding  model  development  is  that  short- 
period  breaking  waves  during  storm  events  are  the  dominant  cause  of  beach 
change  In  the  surf  zone  and  thus  of  bar  formation  and  movement. 

Lone- period  waves 

29.  Investigations  on  several  coasts  (e.g.,  Guza  and  Thornton  1982, 
Creenwood  and  Sherman  1984,  Wright  and  Short  1984)  have  found  that  significant 
amounts  of  energy  can  be  contained  in  lower  frequencies  of  the  wave  power 
spectra  in  very  shallow  water  and  in  the  swash.  The  presence  of  lower- 
frequency  motion  cannot  be  predicted  at  present,  but  it  appears  to  be  assoc¬ 
iated  with  storms.  It  has  also  been  shown  that  long  waves  and  periodic 
arrival  of  high  (short-period)  wave  groups  can  be  responsible  for  as  much  or 
more  sediment  suspension  than  individual  short-period  waves  (e.g.,  Downing 
1984;  Sternberg,  Shi,  and  Downing  1984;  Beach  and  Sternberg  1987)  In  shallow 
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portions  of  the  surf  zone.  Thus,  if  long-period  waves  or  wave  groups  are 
present,  their  contribution  to  the  total  sediment  transport  Is  expected  to  be 
significant.  An  equally  if  not  more  significant  role  played  by  long  waves  and 
wave  groups  than  in  direct  movement  of  sediment  is  their  contribution  to 
elevating  the  water  level  near  the  shore.  This  allows  short-period  waves  to 
attack  higher  on  the  beach  than  might  be  expected. 

30.  Quantitative  engineering  relationships  for  predicting  sediment 
transport  produced  by  long  waves  and  groups  have  not  yet  been  developed. 
Therefore,  their  contribution  is  not  directly  represented  in  the  model.  When 
such  transport  relations  become  available,  it  should  be  possible  to  incor¬ 
porate  them  in  SBEACH  if  they  can  be  expressed  in  a  macroscale  framework. 

Field  data 

31.  Field  data  from  the  FRF  was  used  to  verify  the  numerical  model, 
primarily  with  regard  to  bar  movement  (Report  1) .  Temporal  changes  in 
shoreline  position  at  the  FRF  are  relatively  small,  making  the  data  poten¬ 
tially  unsuitable  for  testing  the  predictive  capabilities  of  the  model  for 
evaluating  dune  and  foreshore  change.  In  the  present  report,  dune  and 
foreshore  change  predictions  are  evaluated  using  two  recently  acquired  field 
data  sets  augmenting  field  verification  given  in  Report  1. 

Cross-Shore  Transport  Direction  and  Distribution 

32.  Cross-shore  movement  of  sediment  is  governed  by  properties  of  the 
velocity  field  and  sediment  concentration.  Thus,  a  first-principles  descrip¬ 
tion  of  the  waves,  velocity,  and  concentration  is  necessary  not  only  across 
the  surf  zone  but  also  through  the  water  column  at  every  point  on  the  calcula¬ 
tion  grid.  Presently,  such  a  detailed  prediction  is  not  possible  in  numerical 
models  of  beach  profile  change  for  engineering  use  to  describe  the  evolution 
of  natural  beaches.  However,  Dally  and  Dean  (1984)  and,  more  recently, 
Roelvlnk  and  Stive  (1989)  have  had  success  with  first-principles  type  models 
in  reproducing  specific  events  of  profile  change. 

33.  Sediment  concentration  in  the  surf  zone  is  closely  related  to  the 
generation  of  turbulent  motion,  which  depends  on  wave  breaking.  Accordingly, 
it  is  plausible  that  sediment  concentration,  and  thus  the  amount  of  material 
available  for  transport,  is  closely  related  to  wave  energy  dissipation 
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produced  by  wave  breaking  (Dean  1977,  Krlebel  1982,  Moore  1982,  Krlebel  and 
Dean  1985).  This  macroscale  approach  Is  used  In  SBEACH.  In  the  model, 
prediction  of  the  cross-shore  velocity  field  is  avoided  by  expressing  the 
cross-shore  sand  transport  rate  directly  in  terms  of  wave  and  profile  proper¬ 
ties.  The  direction  of  transport  is  predicted  from  an  empirical  criterion 
derived  from  the  large  wave  tank  data  and  verified  with  field  data,  and  the 
magnitude  of  transport  in  the  main  portion  of  the  surf  zone  is  related  to  the 
wave  energy  dissipation  per  unit  water  volume.  This  simplified  approach 
implies  a  transport  race  distribution  directed  either  onshore  or  offshore 
along  the  entire  profile  at  a  specific  instant  In  time. 

Cross-shore  transport  direction 

34.  Criteria  for  predicting  whether  a  beach  will  erode  or  accrete 
through  cross-shore  sand  transport  processes  have  been  suggested  by  a  number  of 
authors  (see  Larson  and  Kraus  1989a) .  Most  criteria  refer  to  general  types  of 
profile  configuration,  such  as  winter/summer,  storm/normal,  bar/berm,  dlsslpa- 
tive/reflectlve ,  which  evolve  under  certain  wave  and  sediment  conditions, 
whereas  some  authors  have  focused  on  predicting  the  transport  direction.  In 
general,  there  is  no  contradiction  between  the  two  types  of  criteria,  and  they 
can  be  regarded  as  equivalent  in  a  macroscopic  sense.  A  profile  on  which  bars 
develop  has  transport  directed  predominantly  offshore,  whereas  a  profile 
exhibiting  berm  build-up  mainly  experiences  onshore  sand  movement.  However,  If 
a  criterion  for  predicting  transport  direction  is  expressed  In  terms  of  local 
wave  properties,  there  could  be  a  significant  difference  between  transport 
direction  and  the  overall  profile  type  that  evolves. 

35.  Criteria  for  predicting  profile  type  or  overall  transport  direction 
typically  include  one  parameter  for  characterizing  the  incident  wave  condition 
and  another  parameter  Involving  some  property  of  the  sediment  (grain  size  or 
fall  speed)  and/or  beach  slope.  Deepwater  wave  steepness  Hq/Lq  Is  the  most 
commonly  used  descriptor  of  the  waves  (Waters  1939,  Rector  1954,  Iwagakl  and 
Node  1962,  Nayak  1970,  Dean  1973,  Sunamura  and  Horlkaws  1974),  whereas  several 
different  parameters  for  describing  sediment  properties  have  been  proposed, 
such  as  DS0Ao  (Rector  1954),  H0/DJO  (Iwagakl  and  Noda  1962),  and  «w/gT 
(Dean  1973) ,  In  which 
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H0  -  deepwater  wave  height,  m 
L0  ~  deepwater  wavelength,  a 
Di0  «  median  grain  size,  mm 
w  -  sand  fall  speed,  m/sec 
g  «■  acceleration  due  to  gravity,  m/sec2 
T  -  wave  period,  sec 

In  development  of  SBEACH,  the  deepwater  wave  steepness  Hq/Lq  and  the  para¬ 
meter  Hc/wT  ,  called  the  dimensionless  fall  speed  or  fall  speed  parameter, 
were  found  to  give  a  good  distinction  between  profiles  exhibiting  mainly  bar 
and  berm  formation.  The  fall  speed  parameter  also  emerged  naturally  in  the 
course  of  regression  analysis  relating  morphodynamic  properties  of  the  profile 
to  the  incident  waves  (Larson  1988;  Larson,  Kraus,  and  Sunamura  1988;  Larson 
and  Kraus  1989a) . 

36.  Figure  1  shows  the  criterion  for  distinguishing  erosion  and 
accretion  developed  for  data  from  the  large  wave  tank  experiments.  Separation 
of  eroded  and  accreted  (bar  and  berm)  profiles  is  achieved  by  the  equation 
expressing  the  solid  line  drawn  in  the  figure  and  given  by 


In  which  M  -  0.00070  is  an  empirically  determined  coefficient.  If  the  left 
side  of  Equation  1  is  less  (greater)  than  the  right  side,  the  profile  Is 
predicted  to  erode  (accrete).  Wave  steepness  describes  the  wave  asymmetry, 
whereas  wave  height  and  period  appearing  In  the  fall  speed  parameter  account 
for  the  absolute  magnitudes  of  those  quantities.  The  fall  speed  accounts  for 
the  settling  characteristics  of  sand  particles.  Seville  (19S7)  showed  on  the 
basis  of  the  early  CE  experiments  that  the  magnitude  of  the  wave  height 
controls  erosion  and  accretion  in  addition  to  the  wave  steepness. 

37.  Viewing  Figure  1,  it  la  seen  that  higher  waves  and  shorter  periods 
tend  to  promote  erosion.  Coarser-grained  beaches,  for  which  the  fall  speed  Is 
greater,  tend  to  more  readily  accrete  than  do  finer-grained  beaches  exposed  to 
the  same  wave  conditions.  These  predictions  conform  with  field  observations. 
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Figure  1.  Criterion  for  predicting  beech  erosion  end  accretion 
baeed  on  monochromatic  laboratory  wave  data 


38.  In  describing  profile  change  In  the  wave  tank  studies,  specifics* 
tlon  of  wave  and  sediment  properties  to  be  used  In  Equation  1  Is  unambiguous 
because  of  the  fixed  wave  height  and  period  and  well-defined  sand  size. 
However,  In  order  to  apply  the  equation  to  a  natural  beach,  appropriate  wave 
and  sediment  characteristics  must  be  available.  Field  data  were  collected 
from  various  Investigations  around  the  world  to  establish  the  proper  statis¬ 
tical  descriptors  to  be  used  In  Equation  1.  The  data  used  are  summarized  in 
Sunamura  (1980a)  and  Seymour  (1986),  where  movement  of  the  deeper  contours  was 
used  to  evaluate  profile  change  and  transport  direction  In  the  latter  investi¬ 
gation.  Report  1  gives  a  more  complete  description  of  the  field  data  used  in 
evaluating  the  performance  of  Equation  1. 

39.  Three  deepwater  statistical  vave  heights  (mean  wave  height  R0  , 

root-mean-square  wave  height  ,  significant  wave  height  Hso)  were  tested 

In  Equation  1  to  obtain  the  best  division  between  accretionary  and  erosional 
profiles.  Although  reasonable  delineation  could  be  obtained  with  each  of  the 
wave  heights  by  adjustment  of  the  value  of  the  empirical  coefficient  M  ,  use 
of  mean  wave  height  provided  the  best  delineation  and  gave  the  same  value  of 
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M  as  found  for  the  large  wave  tank  data.  Figure  2  displays  In  analogy  to 
Figure  1  the  criterion  (solid  line)  together  with  the  data  from  natural 
beaches.  If  significant  wave  height  was  reported  In  the  field  data,  the  aean 
field  data,  the  mean  wave  height  was  calculated  as  R0  -  0.626  Hso  assuming  a 
Rayleigh  distribution.  The  aean  or  peak  period  of  the  spectra  was  chosen  as 
the  representative  period  to  be  used  In  Equation  1  together  with  median  grain 
size  of  the  beach  sedinent . 

40.  Mlaura,  Otsuka,  and  Vatanabe  (1986)  compared  profile  change  versus 
transport  direction  and  rate  produced  In  a  snail  wave  tank  In  comparable  cases 
using  regular  and  Irregular  waves.  They  observed  that  the  criterion  developed 
by  Sunanura  and  Horlkava  (1974)  gave  the  best  delineation  between  erosion  and 
accretion  If  aean  wave  height  was  used  In  analyzing  the  Irregular  wave  tests. 
Profile  change  proceeded  at  a  slower  rate  for  Irregular  waves,  and  this  was 
attributed  to  the  presence  of  both  accretionary  and  eroslonal  wave  components 
In  a  wave  train.  Their  observation  supports  the  nethod  of  super lnposlng  the 
effects  of  Individual  waves  to  obtain  profile  evolution  for  irregular  waves  In 
SBEACH .  Kraus  and  Horlkava  (1990)  further  discuss  the  difference  In  sand 
transport  under  regular  and  Irregular  waves. 
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Cross-shore  sand  transport  rate  distribution 

41.  Cross -shore  sand  transport  rate  distributions  were  calculated  from 
consecutive  profiles  In  time  by  Integrating  the  mass  conservation  equation. 
This  technique  has  been  applied  by  van  Hijum  (1974,  1976),  Hattori  and 
Kawamata  (1980),  and  Shimizu  et  al.  (1985).  The  transport  rate  obtained 
constitutes  an  average  net  rate  for  the  time  period  between  the  tvo  profile 
surveys.  Report  1  describes  a  wide  range  of  analyses  of  cross -shore  transport 
race  characteristics  and  transport  rate  formulas  derived  on  the  basis  of 
physical  considerations  and  observations  from  the  data.  In  the  following, 
some  Important  aspects  of  the  cross -shore  sand  transport  rate  are  discussed 
regarding  typical  distribution  shapes  in  time  and  space. 

42.  An  important  property  of  the  cross -shore  sand  transport  rate  as 
calculated  from  measured  profile  change  is  the  temporal  decay  of  the  rate  of 
change  for  profiles  exposed  to  constant  waves  and  water  level.  The  calculated 
maximum  sand  transport  rate  was  more  chan  one  order  of  magnitude  smaller  at 
the  end  of  a  run  as  compared  to  the  earliest  part  of  the  run.  A  small  amount 
of  transport  activity  still  occurred  at  the  end  of  the  runs,  attributed  In 
part  to  profile  change  Introduced  by  small  variations  in  forcing  conditions 
and  randomness  in  fluid  and  sand  particle  motion.  Thus,  In  fitting  empirical 
functions  by  least-squares  to  calculated  maximum  transport  rates,  an  inverse 
dependence  on  time  gave  a  better  description  than  an  exponential  decay  because 
of  the  finite  amount  of  transport  prevailing  at  the  end  of  the  run. 

43.  By  calculating  the  cross -shore  transport  rate  distribution  using 
Initial  and  final  profiles  of  a  run,  termed  equilibrium  distributions,  a 
picture  of  the  overall  profile  response  is  obtained.  Equilibrium  transport 
rate  distributions  were  classified  as  three  main  types  (Larson  and  Kraus 
1989a) .  Figure  3  Illustrates  che  characteristics  of  these  types  by  displaying 
the  calculated  distributions  from  selected  large  wave  tank  cases.  The  equili¬ 
brium  distribution  types  may  be  summarized  as: 

g.  Eroslonal  (Type  E);  transport  directed  offshore  along  the 
entire  profile. 

fa.  Accretionary  (Type  A);  transport  directed  onshore  along  the 
entire  profile. 

£.  Nixed  Accreclonary  and  Eroslonal  (Type  AE);  transport  directed 
onshore  along  the  seaward  part  of  the  profile  and  transport 
directed  offshore  along  the  shoreward  part  of  the  profile. 
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44.  A  majority  of  cases  fell  In  the  classifications  of  Type  E  or  Type 
A,  indicating  predominantly  unidirectional  transport  during  the  run.  Cases 
having  Type  AE  distributions  correspond  to  profiles  where  little  change  took 
place,  implying  a  transport  direction  which  was  not  so  distinct.  If  sig¬ 
nificant  profile  change  occurred,  transport  was  directed  primarily  onshore  or 
offshore.  Thus,  in  the  numerical  model,  which  has  as  a  primary  objective 
simulation  of  events  producing  significant  profile  change,  it  is  reasonable  to 
assume  a  unidirectional  transport  rate  distribution  as  a  first  approximation. 

45.  Sand  transport  rate  distributions  calculated  from  consecutive 
profile  surveys  did  not  show  as  clear  a  trend  as  did  the  equilibrium  distribu¬ 
tions,  particularly  for  consecutive  surveys  taken  at  later  times.  In  the 
beginning  of  a  run,  when  large  profile  change  took  place,  the  transport  rate 
distribution  was  typically  unidirectional.  Complexity  in  transport  rate 
distributions  at  later  times  probably  reflects  a  close  balance  among  the  many 
factors  contributing  to  sand  movement,  for  which  small  instabilities  away  from 
equilibrium  could  produce  relatively  large  localized  contributions  to  the 
transport.  In  order  to  model  beach  change  in  situations  where  the  beach  is  In 
near-equilibrium  with  the  water  level  and  incident  waves,  prediction  of  more 
complex  distributions  will  probably  have  to  be  Included. 

Tr»nfp<?n_Fegl<?ns,  snflJTrangjprt  Rates 

46.  Report  1  describes  analyses  of  individual  net  sand  transport  rate 
distributions  with  the  aim  of  relating  the  local  transport  race  to  wave  and 
sediment  characteristics.  Research  in  nearshore  wave  dynamics  has  identified 
regions  of  different  hydrodynamic  properties  associated  with  wave  shoaling, 
wave  decay,  and  intensity  of  small-  and  large-scale  motion.  Svendsen,  Madsen, 
and  Buhr  Hansen  (1978)  defined  an  outer  and  inner  region  of  the  surf  zone 
according  to  the  character  of  wave  and  wave -induced  water  motion.  The  outer 
region,  located  immediately  shoreward  of  the  wave  break  point,  is  distin¬ 
guished  by  a  rapid  transition  in  wave  height  with  large-scale  flow  patterns. 

At  some  location  landward  of  the  break  point,  the  large-scale  flow  decays  and 
turns  into  small-scale  turbulent  fluctuations;  simultaneously  the  wave  changes 
character  and  becomes  bore-like.  The  region  shoreward  of  this  point  is  known 
as  the  inner  region  and  extends  to  the  location  where  wave  runup  begins. 
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47.  The  concept  of  two  different  regions  in  the  surf  zone  had  already 
been  suggested  by  Miller  (1976),  who  distinguished  a  vortex  and  a  bore  region 
(see  also  Basco  1985  and  Jansen  1986).  Miller  also  noted  that  the  main 
difference  between  plunging  and  spilling  breakers  was  the  generation  of 
vortices  in  the  outer  region.  In  the  inner  (bore)  region,  flow  conditions 
were  similar  and  Independent  of  breaker  type,  as  also  noted  by  Svendsen, 
Madsen,  and  Buhr  Hansen  (1978).  Svendsen  (1987)  analyzed  flow  measurements  to 
determine  the  turbulent  kinetic  energy  in  the  surf  zone.  He  found  that  the 
variation  over  depth  of  the  time-mean  turbulent  kinetic  energy  was  small 
(Peregrine  and  Svendsen  1978),  and  he  attributed  this  to  a  strong  mixing 
created  by  large-scale  vortices.  Skjelbreia  (1987)  studied  breaking  of 
solitary  waves  in  the  laboratory  and  defined  four  different  zones  of  wave 
transformation:  gradual  and  rapid  shoaling,  and  gradual  and  rapid  decay. 

48.  The  identification  of  regions  with  different  wave  characteristics 
in  the  nearshore  stimulated  a  separate  analysis  of  the  calculated  net  trans¬ 
port  rates  in  these  regions.  An  Important  assumption  underlying  the  numerical 
model  SBEACH  is  the  significance  of  breaking  and  broken  waves  in  determining 
the  cross-shore  sand  transport  rate.  Quite  different  transport  character¬ 
istics  should  be  found  inside  and  outside  the  surf  zone,  making  it  reasonable 
to  assume  that  the  transport  relationships  are  different  in  a  macroscopic 
sense.  Furthermore,  the  dynamics  of  transport  in  the  swash  are  governed  by 
properties  of  the  uprush  bore  and  the  backwash,  sediment  characteristics 
(grain  size,  porosity,  grain  density,  etc.),  and  beach  slope,  which  makes  it 
logical  to  define  the  swash  zone  as  another  distinct  region  of  transport. 

49.  Baaed  on  the  division  of  the  profile  applied  in  nearshore  wave 
dynamics  and  the  physical  characteristics  of  sediment  transport  under  various 
flow  conditions,  four  different  zones  of  transport  (Figure  4)  were  Introduced 
(Larson,  Kraus,  and  Sunamura  1968;  Larson  and  Kraus  1989a).  These  zones  are: 

g.  Zone  I:  From  the  seaward  depth  of  effective  sand  transport 

to  the  break  point  (pre-breaking  zone). 

fe.  Zone  II:  From  the  break  point  to  the  plunge  point  (breaker 
transition  zone). 

£.  Zone  III:  From  the  plunge  point  to  the  point  of  wave  reforma¬ 
tion  or  to  the  swash  zone  (broken  wave  zone). 

&.  Zone  IV:  From  the  shoreward  boundary  of  the  surf  zone  to  the 

shoreward  limit  of  runup  (swash  zone). 
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PLUNGE  BREAK 
POINT  POINT 


SWASH  BREAKER  TRANSITION 

ZONE  BROKEN  WAVE  ZONE _ ^  ZONE  ^ _ PREBREAKING  ZONE 


Figure  4,  Principal  zones  of  cross-shore  transport 

50.  The  zone  located  between  che  break  point  and  the  plunge  point  is 
analogous  to  the  outer  region  (vortex  region)  proposed  by  Svendsen,  Madsen, 
end  Buhr  Hansen  (1978).  A  certain  distance  is  required  after  breaking  before 
turbulent  conditions  are  approximately  uniform  through  che  water  column. 
Assuming  the  cross-shore  transport  rate  is  proportional  to  energy  dissipation 
per  unit  water  volume,  the  recognition  of  this  zone  is  important  and  physi¬ 
cally  motivated.  In  this  report  only  empirically-based  transport  relation¬ 
ships  for  the  zones  derived  from  the  large  wave  tank  data  will  be  presented. 
These  formulas  are  used  in  the  numerical  model  to  determine  the  transport  rate 
distribution  across-shore  and  are  discussed  in  detail  in  Report  1. 

51.  The  transport  rate  relationships  derived  based  on  physical  con¬ 
siderations  and  analysis  of  the  large  wave  tank  data  may  be  summarized  for  the 
four  different  zones  as 


Zone  I : 

q  -  Rb 

e-A,(x-xb) 

Xb  c  X 

(2) 

Zone  II: 

q  -  qp 

e-X2(x-xp) 

xp  <  x  <;  xb 

(3) 
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Zone  III:  q  - 


0  >[°eq 
0  S[°eq 


t  dhl 
K  dxj 

t 

K  dxj 


(<*) 


Zone  IV: 


x„  I  x  i  x. 

es  ~P  (5) 

Xj.  <  X  <  xeg 


where 

q  -  net  cross -shore  sand  transport  rate,  mJ/n/sec 
x  -  cross-shore  coordinate  directed  positive  offshore,  d 
X,2  -  spatial  decay  coefficients  in  Zones  I  and  II,  1/m 
K  -  sand  transport  rate  coefficient,  m*/N 

D  -  wave  energy  dissipation  per  unit  water  volume,  Nm/ms/sec 
Deq  -  equilibrium  wave  energy  dissipation  per  unit  water  volume, 

Nm/m3/sec 

t  -  slope-related  sand  transport  rate  coefficient,  m2/sec 
h  -  water  depth,  m 

The  subscripts  b  ,  p  ,  es  ,  and  r  stand  for  quantities  evaluated  at  the 
break  point,  plunge  point,  end  of  the  surf  zone,  and  runup  limit,  respec¬ 
tively.  Different  spatial  decay  coefficients  are  used  In  Zones  I  and  II, 
denoted  by  the  subscripts  1  and  2 ,  to  describe  the  decrease  in  sand  transport 
rate  with  distance. 

52.  Empirical  expressions  for  the  spatial  decay  coefficients  were 
derived  from  large  wave  tank  data  (Report  1}.  The  decay  coefficient  in  Zone  I 
was  empirically  related  to  median  grain  size  and  breaking  wave  height  as 


Xj  - 


(6) 


For  Zone  II,  limited  data  suggested  that 

X2  -  0.2  Xj  (7) 
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In  calculating  the  transport  rate  in  Zones  I  and  II,  the  transport  rate  is 
first  determined  at  the  plunge  point  from  Equation  4,  and  then  the  exponential 
decay  rates  are  applied  seaward  in  the  respective  zones. 

53.  The  cross -shore  transport  rate  in  zones  of  fully  broken  waves  is 
related  to  the  wave  energy  dissipation  per  unit  water  volume.  This  type  of 
transport  formula  has  previously  been  used  by  Moore  (1982),  Krlebel  (1982, 
1986),  and  Kriebel  and  Dean  (1985).  Analysis  of  the  large  wave  tank  data 
substantiated  this  type  of  transport  relationship  in  zones  of  broken  waves 
(Larson  1988,  Larson  and  Kraus  1989a).  The  transport  relationships  in  the 
other  zones  are  empirical  and  based  directly  on  the  data  from  the  wave  tank 
experiments.  For  example,  in  the  swash  zone  the  transport  rate  simply 
decreases  linearly  from  the  end  of  the  surf  zone  to  the  runup  limit.  This 
property  of  the  transport  rate  distribution  appeared  for  many  large  wave  tank 
cases,  occurring  with  both  onshore  and  offshore  transport.  Similar  character¬ 
istics  have  been  noted  on  a  natural  beach  where  the  foreshore  changed 
uniformly  during  erosional  and  accretionary  events  (Seymour  1987).  However, 
no  direct  measurements  have  been  made  of  transport  rates  on  the  foreshore,  and 
this  region  of  the  profile  deserves  further  study. 
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PART  III:  WAVE  MODEL 


54.  The  wave  model  provides  Input  to  calculate  cross -shore  sand  trans¬ 
port  rates  from  which  profile  change  is  obtained.  In  this  chapter,  an 
overview  is  given  of  the  capabilities  and  limitations  of  the  wave  model 
including  a  summary  of  the  governing  equations.  The  equations  are  expressed 
in  an  explicit  finite-difference  form  and  presented  together  with  the  numer¬ 
ical  solution  method.  Finally,  the  interaction  of  the  wave  model  with  other 
modules  in  the  full  profile  change  model  is  discussed  with  reference  to  input 
and  output  requirements. 

55.  The  wave  model  is  a  generalized  version  of  the  model  described  in 
Report  1.  Significant  improvements  include: 

g.  Use  of  complete  linear  wave  theory  expressions  rather  than 
shallow  water  wave  approximations . 

fe.  Addition  of  an  option  for  specifying  an  arbitrary  incident  wave 
angle  with  respect  to  the  bottom  contours.  The  bottom  contours 
are  assumed  to  be  straight  and  parallel. 

£.  Use  of  a  more  general  numerical  solution  scheme  based  on  the 
energy  flux  equation. 

Model  Capabilities  and  Limitations 


56.  Linear  wave  theory  is  used  to  determine  wave  characteristics 
across-shore  from  deep  water  or  a  specified  water  depth  offshore  to  the  break 
point.  Shoreward  of  the  break  point,  a  slightly  generalized  form  of  the  wave 
decay  model  of  Dally  (1980)  and  Dally,  Dean,  and  Dalrymple  (1984)  calculates 
the  wave  height  distribution.  The  point  of  incipient  wave  breaking  is 
determined  from  an  empirical  criterion  expressed  in  terms  of  the  surf  simi¬ 
larity  parameter  defined  with  the  deepwater  wave  steepness  and  the  local  beach 
slope  seaward  of  the  breakpoint,  evaluated  over  one-third  of  the  wavelength. 

57.  An  input  wave  height,  period,  and  Incident  wave  angle  must  be  given 
either  in  deep  water  or  at  a  specified  depth.  Vave  information  is  transformed 
by  shoaling  and  refraction  over  straight  and  parallel  bottom  contours  to  the 
most  seaward  point  of  the  model  grid.  Energy  dissipation  in  the  viscous  layer 
along  the  bottom  is  neglected  since  the  effect  on  wave  height  is  negligible  in 
comparison  to  the  amount  of  energy  dissipation  in  the  surf  zone  (Dally  1980). 


58.  Calculation  of  the  wave  height  across -shore  starts  In  the  moat 
seaward  calculation  cell  and  proceeds  shoreward  cell  by  cell.  The  local  wave 
height  and  angle  with  respect  to  the  bottom  contours  are  determined  from 
conservation  of  energy  flux  and  Snell’s  Law.  In  regions  of  wave  breaking, 
energy  dissipation  due  to  Increased  turbulence  In  the  water  column  Is  Included 
In  the  energy  flux  equation.  It  Is  possible  to  specify  bottom  contours  with 
an  orientation  that  Is  not  parallel  to  the  baseline  (y-axls). 

59.  Energy  dissipation  In  the  surf  zone  Is  assumed  proportional  to  the 
excess  energy  flux  beyond  a  stable  energy  flux  which  Is  a  function  of  water 
depth  (Horlkawa  and  Kuo  1966) .  Wave  properties  In  the  surf  zone  are  calcu¬ 
lated  using  linear  wave  theory  without  Invoking  shallow  water  approximations. 
The  wave  height  calculation  Is  carried  out  to  a  water  depth  which  is  less  than 
the  depth  where  the  surf  zone  ends  in  accordance  with  the  definition  of  the 
different  hydrodynamic  regions  in  the  nearshore.  The  depth  at  the  shoreward 
end  of  the  surf  zone  Is  roughly  equal  to  the  seaward  limit  of  the  backrush. 

60.  The  wave  model  allows  wave  reformation  to  take  place  and  the  possi¬ 
bility  of  multiple  break  points  across  the  nearshore  zone.  Wave  reformation 
occurs  If  the  local  wave  height  reaches  the  stable  wave  height.  Implying  no 
wave  energy  dissipation.  After  energy  dissipation  ceases,  shoaling  and 
refraction  become  dominant  until  wave  breaking  occurs  and  dissipation  Is 
initiated  once  more. 

61.  Generalizing  the  breaker  decay  model  of  Dally  (1980)  and  Dally, 
Dean,  and  Dalrymple  (1984),  the  two-dimensional  equation  for  conservation  of 
energy  flux  Incorporating  energy  dissipation  associated  with  wave  breaking  may 
be  written 

4  (F  cose)  +  4  (F  sine)  -  £  (F  -  Fs)  (8) 

Sx  3y  h 

where 

F  -  wave  energy  flux,  Nm/m/sec 

8  •  wave  angle  with  respect  to  the  bottom  contours 
x  —  cross-shore  coordinate  positive  In  the  seaward  direction,  m 
y  -  longshore  coordinate,  m 
x  -  wave  decay  coefficient 
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h  -  water  depth,  m 

Fs  -  stable  wave  energy  flux,  Nm/m/sec) 
The  wave  energy  flux  is  given  by 


F  -  ECS 

where 

E  -  wave  energy  density,  Nm/m4 
Cg  -  wave  group  speed,  a/s 

The  wave  energy  density  is  written  using  linear  wave  theory  as 

E  -  i  peh4 


(9) 


(10) 


where 

H  -  wave  height,  m 

p  -  density  of  water,  Vg/m* 

g  -  acceleration  due  to  gravity,  m/sec* 

Assuming  the  wave  conditions  to  be  uniform  alongshore  and  the  bottom  contours 
to  be  straight  and  parallel.  Equation  8  may  be  reduced  to 


3  x  _  , 
—  (F  cos6)  -  -  (F  -  Fs) 
3x  h 


(ID 


62.  Wave  group  speed  Is  related  to  wave  phase  speed  C  through  a 
factor  n  which  Is  a  function  of  the  water  depth  and  the  wavelength  L  (or 
wave  period  T  )  as 


vhere 


n 


1 

2 


1  + 


(13) 
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The  wave  phase  speed  is  determined  through  the  dispersion  relationship 


C 


Co 


where 


(U) 


(15) 


in  which  Cc  is  the  wave  phase  speed  in  deep  water.  Equation  14  is  an 
implicit  equation  which  must  be  solved  numerically  or  by  rational  approxima¬ 
tions  since  C  contains  a  dependence  on  the  wavelength,  that  is,  C  —  L/T  . 

63.  The  wave  decay  coefficient  controls  the  rate  of  energy  dissipation 
whereas  the  stable  energy  flux  determines  the  amount  of  energy  dissipation 
necessary  for  stable  conditions  to  occur  once  breaking  is  initiated.  In  this 
respect,  the  stable  wave  condition  refers  to  a  state  in  which  energy  dissipa¬ 
tion  during  breaking  ceases,  allowing  waves  to  reform.  The  stable  energy  flux 
may  be  expressed  as 


Fs  -  Escg  U«> 

where  Es  -  stable  wave  energy  density,  Nm/m2.  The  stable  wave  energy  flux 
corresponds  to  a  stable  wave  height  which  is  a  function  of  the  water  depth 

Hs  -  Th  (17) 

where  T  -  stable  wave  height  coefficient. 

64.  Two  empirical  coefficients  enter  in  the  breaker  decay  model,  k 
and  T  .  However,  the  values  of  these  coefficients  seem  to  vary  little 
over  a  wide  range  of  conditions,  both  in  the  laboratory  and  on  natural 
beaches.  On  the  basis  of  small  and  large-scale  tank  data,  Dally,  Dean,  and 
Dalrymple  (1984)  recommended  the  values  of  k  -  0.15  and  T  -  0.40  for  use 
as  averages  for  a  varying  bottom  slope.  Ebersole  (1987)  evaluated  the 
performance  of  the  breaker  decay  model  using  field  measurements  and  found  good 
agreement  for  engineering  applications. 
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65.  Report  1  verified  the  applicability  of  the  breaker  decay  model  for 
a  barred  beach  profile  using  the  large-wave  tank  measurements  of  Kajima  et  al. 
(1983).  Values  obtained  for  the  two  empirical  coefficients  were  In  general 
agreement  with  values  obtained  by  the  aforementioned  authors.  An  apparent 
dependence  on  beach  slope  for  the  coefficients  was  noted  but  no  specific 
relationships  could  be  developed. 

66.  Kraus  and  Larson  (In  preparation)  evaluated  the  performance  of  the 
breaker  decay  model  for  waves  approaching  the  beach  at  an  angle  with  respect 
to  the  bottom  contours  both  for  laboratory  (Vlsser  1982)  and  field  data  (Vu, 
Thornton,  and  Cuza  1985) .  In  these  cases  the  coefficients  were  fixed  at  the 
standard  values  of  k  -  0.15  and  T  -  0.60.  in  all  cases,  the  breaker  decay 
model  produced  satisfactory  agreement.  For  the  field  data,  a  statistical 
approach  was  used  In  which  wave  heights  were  randomly  selected  from  a  Rayleigh 
distribution  for  a  large  number  of  waves,  and  the  calculated  wave  height 
distributions  were  averaged  to  obtain  a  representative  distribution. 

67.  As  waves  propagate  toward  shore,  a  flux  of  momentum  (radiation 
stress)  arises  which  causes  displacement  of  the  mean  water  elevation  if 
changes  In  this  flux  occur  due  to  shoaling  or  breaking.  Seaward  of  the  break 
point  shoaling  produces  an  Increase  In  wave  height  which  causes  a  correspond¬ 
ing  Increase  In  the  momentum  flux.  This  flux  increase  is  balanced  by  lowering 
of  the  mean  water  elevation,  called  setdown  (Longuet-Hlggins  and  Stewart 
1963).  Inside  the  surf  zone,  as  waves  continue  to  break  and  decrease  In 
height,  the  momentum  flux  decreases  and  an  increase  In  mean  water  elevation 
occurs  (setup). 

68.  Displacement  of  the  mean  water  surface  (setup  or  setdown)  may  be 
determined  from  the  momentum  equation 


£fxx 

ax 


3t, 

pg(h  +  ■n)  — 
3x 


(18) 


where 

Sgg  -  radiation  stress  component  directed  onshore,  N/m 
q  -  wave  setup  or  setdown,  m 

The  radiation  stress  component  Sxx  Is  given  by  linear  wave  theory  for  an 
arbitrary  wave  angle  of  Incidence  as 
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Sxx  -  \  PS”2  [n(c0sS6  +  X>  *  i] 


(19) 


69.  The  energy  flux  conservation  equation  (Equation  8),  the  wave 
dispersion  relationship  (Equation  14),  and  the  cross-shore  momentum  equation 
(Equation  18)  are  not  sufficient  to  specify  the  wave  height  across-shore  If 
the  waves  arrive  at  an  oblique  angle  Co  the  bottom  contours.  In  this  case  an 
additional  equation  must  be  solved  which  expresses  the  conservation  of  wave 
number.  The  wave  number  conservation  equation  Is  written  In  the  general  form, 
not  taking  into  account  diffraction,  as 


a 

(sine-) 

a 

(A 

O 

O 

ax 

l  l  J 

3y 

l  Lj 

70.  For  straight  and  parallel  bottom  contours,  Implying  no  alongshore 
variation  in  wave  properties,  Equation  20  may  be  simplified  to  yield 


(21) 


Integration  of  Equation  21  gives  Snell's  law,  which  states  that  the  quantity 
sin8/L  is  constant  across-shore.  Equations  8,  14,  19,  and  21  allow  the  wave 
height  to  be  determined  once  boundary  conditions  have  been  established. 

71.  As  mentioned  previously  the  wave  height,  period,  and  Incident  angle 
are  calculated  at  the  seaward  end  of  the  profile  from  deepwater  or  equivalent 
conditions,  providing  the  necessary  boundary  conditions.  The  end  of  the  grid 
should  be  located  seavard  of  the  break  point  for  all  simulations.  From  the 
wave  properties  it  Is  possible  to  determine  analytically  the  setdown  at  the 
most  seaward  point  (Longuet-Hlgglns  and  Stewart  1962) . 


*1 


xHl 


4L  sinh 


(22) 
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72.  Also,  a  condition  for  the  Initiation  of  breaking  is  needed, 
Indicating  the  location  at  which  energy  dissipation  due  to  breaking  begins  and 
the  wave  decay  coefficient  becomes  different  from  zero.  Breaking  occurs  when 
the  ratio  between  the  wave  height  and  the  water  depth  exceeds  a  certain  value 
which  Is  a  function  of  the  surf  similarity  parameter  expressed  in  terms  of  the 
deepwater  wave  steepness  and  profile  slope  prior  to  breaking. 

T  -  1.14  €»■»  (23) 


where 

t  -  breaker  ratio  -  Hjj/hj, 
b  -  subscript  denoting  wave  breaking 
5  -  surf  similarity  parameter 
The  surf  similarity  parameter  Is  given  by 

i -0.5 

C  -  tang  gM  (24) 

where  tang  »  local  beach  slope  seaward  of  the  break  point. 

73.  Equation  23  was  derived  in  Report  1  from  data  given  by  Kajima  et 
al.  (1983),  but  very  similar  relationships  were  previously  reported  in 
Slngamsettl  and  Wind  (1980)  and  Sunamura  (1980b).  Report  1  used  data  from 
prototype -scale  experiments  whereas  the  two  latter  studies  were  based  on 
small-scale  laboratory  experiments.  Thus,  in  Equation  23  scale  effects  are 
absent  although  it  should  be  observed  that  only  monochromatic  waves  vers  used 
in  the  tests. 

74.  The  laboratory  data  of  Kajima  et  al.  (1983)  Included  cases  with  a 
maximum  wave  period  of  12  sec  although  most  eases  had  periods  In  the  range  of 
3  to  9  sec.  The  limited  data  available  for  longer  period  waves  indicate  that 
Equation  23  should  be  used  with  care.  Also,  since  the  relationship  was 
determined  mainly  from  barred  profiles  exposed  to  monochromatic  waves,  the 
profile  slope  seaward  of  the  break  point  was  steeper  than  what  Is  typically 
found  in  the  field.  Thus,  Equation  23  tends  to  underestimate  the  breaker 
ratio  for  gentle  slopes  in  field  applications. 
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75.  An  expression  for  energy  dissipation  due  to  wave  breaking  is  needed 
in  the  prediction  of  the  cross -shore  sand  transport  rate.  In  the  process  of 
wave  height  decay  after  breaking,  turbulent  motion  is  induced  in  the  surf 
zone.  The  turbulent  kinetic  energy  is  converted  into  heat  through  transfer  of 
energy  from  larger  to  smaller  eddies.  In  the  macroscale  approach  of  the  beach 
profile  model,  it  is  satisfactory  to  derive  energy  dissipation  from  wave 
height  decay  without  resolving  the  Internal  structure  of  this  process. 

However,  the  difference  between  the  production  of  turbulent  kinetic  energy  and 
the  corresponding  dissipation  with  reference  to  the  time  scale  (Roelvlnk  and 
Stive  1989)  is  observed  as  well  as  the  process  of  energy  reordering  (Svendsen, 
Madsen,  and  Buhr  Hansen  1978) . 

76.  The  wave  energy  dissipation  D  per  unit  water  volume  is 

D  -  ^  r—  (F  cosfl)  (25) 

h  dx 

Using  Equation  11  the  following  relationship  for  D  is  obtained 

D  -  p  (F  -  F,)  (26) 

Once  the  wave  height  distribution  is  calculated,  the  energy  dissipation  per 
unit  water  volume  in  the  surf  zone  is  determined  explicitly  from  Equation  26. 
Outside  the  surf  zone  k  Is  set  to  zero  and  no  energy  dissipation  takes  place 
since  bottom  friction  is  neglected. 

Numerical  Solution  Scheme 

77.  An  explicit  finite-difference  scheme  is  used  to  solve  the  equations 
for  determining  cross-shore  wave  height.  In  the  following  section,  an  outline 
is  given  of  the  numerical  solution  scheme  and  the  difference  equations 
employed.  Figure  5  is  a  definition  of  the  grid  and  Illustrates  where  the 
various  quantities  in  the  difference  equations  are  defined. 

78.  In  the  model,  a  quasi-stationary  approach  is  applied,  implying  that 
a  steady-state  solution  of  the  wave  height  distribution  is  determined  at  every 
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Cine  seep.  Propagation  of  individual  waves  is  not  described,  but  it  is 
assumed  that  the  model  time  step  is  significantly  larger  than  the  wave  period. 

79.  For  irregular  wave  conditions,  the  wave  height  used  at  each  time 
step  may  be  regarded  as  a  statistical  measure  which  represents  the  average 
properties  of  the  breaking  waves  during  the  time  step.  This  wave  height  Is 
different  from  the  wave  height  needed  to  reproduce  the  average  wave  height 
distribution  during  a  time  step,  which  Includes  all  broken  and  unbroken  waves 
at  specific  locations.  The  selected  wave  height  parameter  should  represent 
the  average  location  of  the  surf  zone,  giving  a  clear  distinction  between 
regions  of  broken  and  unbroken  waves.  From  experience  with  field  simulations, 
slgnlficsnt  wave  height  Is  satisfactory  in  this  respect  (Report  1). 

80.  Numerical  calculations  start  at  tha  seaward  end  of  the  grid  and 
proceed  onshore  through  an  explicit  solution  scheme  in  which  quantities  known 
at  a  specific  grid  point  are  used  to  determine  corresponding  quantities  at  the 
next  grid  point.  In  the  difference  equations  presented,  the  index  1  is  used 
to  denote  the  number  of  a  specific  grid  point.  A  staggered  grid  Is  used  with 
different  quantities  taken  at  points  In  the  middle  of  calculation  cells  or  at 
the  boundaries  between  cells  (Figure  5).  The  primary  quantity  In  the  middle 
of  a  cell  is  water  depth,  and  these  grid  points  are  termed  h-points.  At  tha 
boundaries  of  calculation  cells,  the  main  quantity  Is  the  cross -shore  trans¬ 
port  rate,  and  quantities  defined  here  are  termed  q-polnts. 


q  i-l  q  i  q  i+l  q  t+2 


Quantities  evoluoted  ot  q-points: 
Quantities  evaluated  ot  h-points: 


Figure  5.  Definition  sketch  of  numerical  grid 
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81.  As  previously  mentioned,  the  wave  height,  period,  and  angle  must  be 
known  at  the  most  seaward  grid  point  (given  or  calculated  from  some  reference 
depth  offshore)  prior  to  starting  the  calculation.  From  these  wave  proper¬ 
ties,  wave  setdown  Is  determined  from  Equation  22.  Profile  depth  Is  given  in 
the  middle  of  each  calculation  cell  whereas  setdown  and  setup  are  calculated 
at  q-polnts.  Depths  at  the  boundaries  of  the  cells  are  given  by  linear 
interpolation.  Also,  energy  flux  and  radiation  stress  In  the  direction  of  the 
waves  are  determined  from  knowledge  of  the  depth  and  wave  quantities  at  the 
first  q-polnt. 

82.  Moving  from  a  specific  q-polnt  to  the  next,  wave  refraction  la 
first  determined  If  the  Incident  wave  approaches  with  an  angle  to  the  bottom 
contours.  From  Equation  21,  the  angle  between  the  wave  crests  and  the  bottom 
contours  at  the  next  shoreward  q-polnt  Is  given  by 

6,  -  arcslnf- - - ]  (27) 

Um  J 

Note  that  grid  point  numbering  Increases  in  the  seaward  direction  since  the  x- 
axls  points  offshore,  but  the  calculation  proceeds  from  offshore  In  the  shore¬ 
ward  direction.  Water  depth  at  grid  point  1  In  Equation  27  Is  corrected 
with  setdown  (setup)  using  the  value  at  grid  point  1+1  to  avoid  use  of  an 
Implicit  or  Iterative  scheme.  Wavelengths  are  calculated  using  the  dispersion 
relationship  (Equation  14)  solved  by  a  Padd  approximation  (Hunt  1979). 

83.  The  next  step  In  the  calculation  Is  to  determine  the  energy  flux 
and  thus  the  wave  height.  Equation  11  Is  vritten  in  difference  form  as 

F‘  "  tF,+1  (co,®m  '  0-5Ac‘)  +  Ac  F*J  (28) 


where 


A 


Cl 


kAx 

hi  +  ^i«i 


(29) 


and  Ax  -  length  step  (m) .  The  stable  wave  energy  flux  is  determined  from 
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F, 


(30) 


s4  -  -  Pg  [r  <hi  +  ’Vi)] 


2  CS»  *  SlM 


An  average  value  for  the  wave  group  speed  Is  used  since  this  quantity  Is  taken 
at  q-polnts,  and  the  stable  wave  energy  flux  is  evaluated  at  h-polnts. 

Seaward  of  the  break  point,  k  is  set  to  zero,  indicating  that  ACj  Is  also 
zero,  since  no  energy  dissipation  due  to  breaking  occurs. 

64.  After  the  energy  flux  has  been  calculated  at  a  specific  point,  the 
corresponding  wave  height  is  determined  from  Equations  9  and  10  as 


Hi 


(31) 


Using  the  wave  height,  radiation  stress  is  calculated  from  Equation  19  and 
setdown  (setup)  is  given  from  Equation  18,  expressed  in  difference  form  as: 


’ll 


-  Vx  + 


<SXX>n.i  '  (SXX>t 
Pg(ht  +  Hi*i) 


(32) 


85.  At  every  calculation  step  a  check  is  made  to  determine  if  wave 
breaking  has  occurred  according  to  Equation  23.  Once  breaking  is  initiated 
the  wave  decay  coefficient  is  given  the  value  x  -  0.15  and  energy  dissipa¬ 
tion  takes  place.  Also,  if  breaking  has  occurred  a  check  is  made  if  wave 
reformation  is  possible,  that  is,  if  the  wave  energy  flux  is  lower  than  the 
stable  energy  flux.  The  stable  wave  height  coefficient  is  set  to  T  -  0.40, 
and  if  there  is  wave  reformation,  x  is  given  the  value  zero.  An  arbitrary 
number  of  surf  zones  with  intermediate  zones  of  wave  reformation  may  occur. 

86.  The  breaker  ratio  is  determined  from  Equation  23  for  the  most 
seaward  break  point.  If  multiple  breaking  occurs,  a  constant  breaker  ratio 
T  -  1-0  is  applied  at  the  more  shoreward  break  points.  The  shoreward 
boundary  of  the  calculation  grid  is  located  at  a  pre-determlned  water  depth 
approximately  corresponding  to  the  seaward  limit  of  the  backrush  or  the 
location  of  a  seawall. 


37 


87.  Wave  height  calculations  are  not  carried  out  using  the  exact 
topography  of  Che  profile;  instead,  a  smoothed  profile  is  used.  The  reason 
for  this  is  that  the  response  of  waves  to  bottom  topography  takes  place  over  a 
certain  distance,  suggesting  that  shoaling  and  refraction  do  not  occur  instan¬ 
taneously.  For  example,  a  very  small  hump  (a  rock)  on  the  bottom  would  not 
cause  a  wave  to  shoal,  which  the  equations  would  predict  if  no  profile 
smoothing  were  applied.  The  number  of  calculation  cells  over  which  the 
smooching  is  performed  (a  moving  average  scheme)  is  determined  based  on  the 
breaking  wave  height  as  3Hb  .  A  simple  predictive  formula  is  used  to 
estimate  the  breaking  wave  height  at  each  time  step  prior  to  determining  the 
wave  height  distribution 


Hb  -  0.53  H0 


-0.2* 


(33) 
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PART  IV:  NUMERICAL  MODEL  OF  STORM- INDUCED  BEACH  CHANGE 


88.  The  model  SBEACH  consists  of  three  calculetion  modules  which  are 
executed  consecutively  at  each  time  step  In  a  simulation.  The  modules 
calculate  wave  height  across-shore,  net  cross-shore  sand  transport  rate,  and 
profile  change,  respectively.  Wave  height  is  used  to  calculate  the  transport 
rate  from  which  profile  change  is  calculated.  The  bathymetry  is  then  revised 
and  these  modules  are  stepped  through  successively  to  calculate  profile  change 
for  the  simulation  period. 

89.  Part  III  describes  the  theoretical  background  and  numerical 
solution  scheme  of  the  wave  module.  Relationships  used  to  determine  the  net 
cross-shore  sand  transport  rate  are  given  in  Part  II.  In  this  chapter,  an 
overview  of  the  capabilities  and  limitations  of  the  profile  change  model  are 
given  with  a  discussion  of  the  input  data  requirements  and  types  of  output 
produced  by  the  model.  Finally,  the  numerical  solution  scheme  for  calculating 
profile  change,  boundary  conditions,  model  stability,  and  typical  time 
interval  and  length  step  are  described. 

Model  Capabilities  and  Limitations 

90.  A  fundamental  assumption  of  the  numerical  model  is  that  breaking 
waves  are  the  major  cause  of  profile  change,  and  thus  the  major  part  of  cross¬ 
shore  sand  transport  takes  place  in  the  surf  Zone.  The  volume  of  sand 
transported  is  determined  by  the  wave  energy  dissipation  per  unit  water  volume 
in  the  surf  zone  (dependent  upon  wave  conditions  and  profile  shape)  and  by 
empirical  functional  relationships  obtained  from  the  large  wave  tank  data  in 
other  zones  along  the  profile.  The  direction  of  transport  is  given  by  an 
empirical  predictive  criterion  (Equation  1)  and  calculation  of  the  fluid 
velocity  field  in  the  surf  zone  is  not  required. 

91.  Wave  height,  wave  period,  and  water  level  can  be  specified  at  each 
time  step  to  account  for  variable  forcing  conditions.  Veter  level  change 
appears  in  the  computations  as  a  translation  of  the  origin  of  the  vertical 
coordinate  axis,  shifted  back  to  the  original  level  prior  to  entering  the  mass 
conservation  equation.  For  a  random  wave  field,  wave  height  and  period  used 
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in  the  calculations  should  be  representative  for  the  current  time  step  with 
reference  to  the  location  of  the  break  point  and  width  of  the  surf  zone. 

92.  The  numerical  model  Incorporates  the  principle  of  superposition  in 
the  sense  that  random  waves,  as  exist  in  the  field,  are  represented  by  a  time 
series  of  monochromatic  waves.  This  assumption  is  appropriate  with  respect  to 
the  engineering  utility  of  the  model  where  a  macroscopic  quasi -stationary 
approach  is  taken  in  calculating  the  cross -shore  wave  height  distribution. 

93.  For  the  large  wave  tank  data,  choice  of  input  wave  height  and 
period  for  the  computations  is  unambiguous  since  waves  of  constant  height  and 
period  were  used.  On  a  natural  beach  the  forcing  parameters  show  appreciable 
variation  in  time.  By  trial  calculations  using  field  data,  it  was  concluded 
that  significant  wave  height  gives  the  best  result  regarding  the  transport 
rate  and  corresponding  profile  change  as  compared  to  mean  or  root -mean- square 
wave  height  (Report  1).  An  intuitive  explanation  for  using  significant  wave 
height  is  that  energy  dissipation,  which  determines  the  transport  rate, 
Increases  with  the  square  of  the  wave  height  Implying  that  higher  waves  give  a 
relatively  larger  contribution.  Thus,  as  a  representative  wave  height,  a 
larger  measure  should  be  used  rather  than  the  mean  wave  height. 

94.  Mlmura,  Otsuka,  Watanabe  (1986)  conducted  small-scale  laboratory 
experiments  using  Irregular  waves  in  order  to  determine  the  representative 
monochromatic  wave  height  for  describing  various  cross-shore  transport 
phenomena.  They  found  that  the  overall  profile  response  as  eroslonal  or 
accretionary  was  best  described  by  the  mean  wave  height.  This  is  in  agreement 
with  the  results  in  Report  1  as  displayed  by  the  criterion  separating  erosion 
and  accretion  (Equation  1  and  Figures  1  and  2).  The  study  by  Mimura,  Otsuka, 
and  Watanabe  did  not  provide  clear  guidance  concerning  choice  of  a  representa¬ 
tive  wave  for  calculating  the  net  cross -shore  transport  rate,  but  both  the 
mean  and  significant  wave  height  appeared  to  be  applicable.  Their  transport 
relationship  (Watanabe  1982)  is  based  on  the  Shields  parameter,  which  Is 
conceptually  and  functionally  different  from  the  wave  energy  dissipation 
approach  used  in  Report  1. 

95.  Median  grain  diameter  for  the  studied  profile  is  a  key  parameter  in 
calculating  the  net  cross-shore  transport  rate  and  resulting  profile  change. 
This  representative  sand  size  determines  the  equilibrium  energy  dissipation 
which  appears  in  the  transport  rate  formula  in  the  zones  of  broken  waves 
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(Equation  4).  Hoore  (1982)  presented  a  relationship  between  and  Djq 

In  graphical  fora  based  on  field  profiles,  and  Dean  (1987)  re-expressed  the 
relationship  In  terms  of  fall  speed.  As  described  In  Report  1,  the  numerical 
model  gave  the  best  overall  fit  for  large  wave  tank  data  If  0.7S  De q  was  used 
In  the  simulations  rather  than  .  This  proved  satisfactory  for  field  data 

as  well. 

96.  In  the  numerical  model  it  Is  possible,  in  principle,  to  specify  a 
different  grain  size,  and  thus  equilibrium  energy  dissipation,  In  every 
calculation  cell.  However,  since  the  movement  of  individual  grains  Is  not 
tracked,  the  median  grain  diameter  will  remain  the  same  throughout  the 
simulation  In  a  specific  cell,  which  Is  unrealistic  considering  grain  move¬ 
ment.  Some  beach  profiles  exhibit  grain  size  characteristics  with  distinctly 
different  regions  across-shore .  If  for  Instance  the  dune  region  has  a  median 
grain  diameter  different  from  the  subaqueous  part  of  the  profile,  It  may  be 
permissible  to  schematically  represent  these  conditions.  In  this  case,  every 
cell  located  In  the  dune  region  Is  assigned  one  grain  size  whereas  the  cells 
located  seaward  are  given  a  different  grain  size.  Using  this  technique.  It  Is 
Implicitly  assumed  that  as  grains  pass  from  the  dune  region  to  the  seaward 
region,  the  transported  amount  exchanged  between  regions  Is  so  small  that  It 
does  not  alter  the  original  distribution  of  grain  size. 

97.  Development  of  the  numerical  model  was  directed  towards  predicting 
profile  change  during  eroslonal  and  accretionary  events.  Two  main  areas  of 
application  of  the  model  are  estimation  of  profile  erosion  and  shoreline 
retreat  during  storm  or  hurricane  events  and  Initial  adjustment  of  beach  fills 
to  changes  In  water  level  and  waves.  Useful  descriptive  quantities  are  eroded 
volume  above  a  certain  contour  and  retreat  distance  of  that  contour.  The 
program  calculates  these  quantities  for  an  arbitrarily  specified  contour  and 
selected  time  step.  These  applications,  together  with  the  properties  of  the 
large  wave  tank  experiments,  emphasized  model  development  focusing  on 
eroslonal  processes. 
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Profile  Change  Model 


98.  Profile  change  Is 
using  the  net  transport  rate 
S.  If  sand  sources  or  sinks 


calculated  from  the  mass  conservation  equation 
distribution  determined  from  Equations  2  through 
are  absent,  the  mass  conservation  equation  Is 


3h  3q 
3t  3x 


(34) 


The  boundary  conditions  are  zero  sand  transport  rate  at  the  runup  limit  and  at 
the  seaward  depth  signifying  the  onset  of  significant  sand  movement.  The 
seaward  boundary  of  the  calculation  grid  (limiting  depth  for  sand  motion)  Is 
effectively  determined  by  the  exponential  decay  In  Zone  I,  taken  to  be  a  point 
where  the  transport  rate  has  decreased  to  a  small  value. 

99.  As  the  wave  height  distribution  is  calculated  across-shore  for  a 
given  time  step,  the  location  of  the  boundaries  between  the  different  sand 
transport  zones  Is  determined.  The  location  of  the  break  point  Is  given 
directly  from  the  wave  height  calculations,  and  the  plunge  point  Is  computed 
as  3Hb  based  on  findings  by  Galvin  (1969)  and  Svendsen  (1987).  The  landward 
boundary  of  the  surf  zone  Is  defined  as  the  point  of  tangency  between  a 
linearly  sloping  beach  face  and  an  equilibrium  profile  for  a  given  median 
grain  diameter,  whereas  the  runup  limit  is  given  by  a  formula  for  predicting 
the  active  profile  height  developed  from  the  large  wave  tank  data  (Report  1). 
The  expression  for  active  profile  height  Is  similar  to  the  formula  provided  by 
Hunt  (1959),  where  the  runup  height  normalized  with  the  incident  wave  height 
is  a  function  of  the  surf  similarity  parameter. 

100.  The  active  subaerial  profile  height  Is  given  by 


ZR  -  1.47  (#-7»  (35) 

The  surf  similarity  parameter  6  Is  defined  In  Equation  24,  and  the  slope  to 
be  used  In  Equation  35  Is  the  average  beach  slope  in  the  surf  zona.  Origi¬ 
nally  Equation  35  was  based  on  the  Initial  slope,  which  has  little  physical 
meaning,  but  application  of  the  model  to  field  data  showed  that  using  the 
average  surf  zone  slope  gave  satisfactory  results.  The  runup  height  has 
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little  influence  on  bar  size  and  movement,  but  it  does  define  the  shoreward 
limit  for  cross-shore  sand  movement  (Report  1).  The  amount  of  erosion  on  the 
foreshore  is  similar  even  for  significant  variations  in  the  runup  height,  but 
the  shape  of  the  foreshore  is  quite  different.  Transport  at  the  landward  end 
of  the  surf  zone  determines  how  much  material  is  transported  from  or  to  the 
foreshore.  Higher  runup  implies  a  larger  area  from  which  material  can  be 
taken,  and  thus  a  different  resulting  foreshore  shape  compared  to  that  given 
by  a  smaller  runup  height. 

101.  From  knowledge  of  the  locations  defining  the  various  transport 
zones,  the  transport  rate  distribution  can  be  calculated  from  Equations  2 
through  5.  Several  break  points  may  occur  along  a  profile  if  wave  reformation 
takes  place,  leading  to  several  zones  of  Type  II  and  III  transport.  In  order 
to  determine  the  transport  rate  distribution,  sand  transport  is  first  calcu¬ 
lated  in  zones  of  fully  broken  waves  (Zone  III)  according  to  Equation  4, 
written  in  difference  form 

<lt  “  J  [(°l  +  ‘  Deq  +  ^h‘  ‘  h»-i)j  <36> 


The  definition  of  the  calculation  grid  is  given  in  Figure  5.  Energy  dissipa¬ 
tion  per  unit  water  volume  is  taken  at  h-polnts  (in  the  middle  of  a  calcula¬ 
tion  cell)  making  it  necessary  to  average  D  when  calculating  the  transport 
rate. 

102.  Calculation  of  the  transport  rate  in  zones  of  fully  broken  waves 
defines  the  boundaries  of  Zone  III,  from  which  the  transport  rate  may  be 
calculated  in  the  other  zones .  After  these  values  are  determined  at  the 
plunge  point  and  at  the  end  of  the  surf  zone,  Equations  2,  3,  and  5  are 
applied  to  completely  define  the  transport  rate  distribution.  Vhen  discre¬ 
tizing  the  mass  conservation  equation,  transport  rate  distributions  from  two 
time  levels  are  used.  The  superscript  k  is  used  to  denote  a  specific  time 
level,  and  the  difference  equation  is  expressed  as 


•  a 


At 


1  fjgl  •  +  q»«i  -  <} 

2  l  Ax  Ax  J 


(37) 
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where  Ac  Is  the  tine  step  (sec).  Since  the  transport  rate  distribution  is 
determined  from  different  relationships  depending  on  the  zone  of  transport, 
its  spatial  derivative  will  generally  be  discontinuous  at  the  boundaries 
between  the  zones.  To  obtain  a  smoother,  nore  realistic  transport  distribu¬ 
tion,  a  3-point  filter  is  applied  tc  the  calculated  transport  rates  by 

q't  -  0.25  qt.,  +  0.50  q,  +  0.25  qln  (38) 

where  the  prime  (')  denotes  the  smoothed  transport  rate. 

103.  In  some  large  wave  tank  cases  the  local  beach  slope  exceeded  the 
angle  of  repose  due  to  constant  waves  and  water  level.  Monochromatic  waves 
caused  erosion  and  accumulation  to  occur  at  the  same  locations  along  the 
profile  for  a  sufficiently  long  duration  to  produce  excessive  profile  steepen¬ 
ing.  Steepest  profile  slopes  occurred  during  eroslonal  conditions  on  the 
shoreward  side  of  the  bar  just  seaward  of  the  berm  crest.  Local  profile 
slopes  exceeding  the  angle  of  repose  rarely  develop  on  a  natural  beach  since 
wave  height,  wave  period,  and  water  level  are  continuously  varying. 

104.  The  concept  of  avalanching,  as  discussed  by  Allen  (1970),  was 
Incorporated  in  the  numerical  model,  and  a  routine  was  Included  to  account  for 
transport  Induced  by  slope  failure .  Allen  recognized  two  different  limiting 
profile  slopes,  the  angle  of  initial  yield  (angle  of  repose)  and  the  residual 
angle  after  shearing.  If  the  local  slope  exceeds  the  angle  of  initial  yield 
4>ly  ,  material  is  redistributed  along  the  slope  and  a  new  stable  slope  is 
attained,  known  as  the  residual  angle  after  shearing  pra  .  From  analysis  of 
maximum  local  profile  slopes  from  the  large  wave  tank  data,  considerably  lower 
values  were  obtained  of  p^y  and  pra  than  reported  by  Allen.  From  his 
experiments  with  natural  sands  (diameters  ranging  from  0.27  to  3.17  mm),  he 
noted  an  angle  of  about  48  deg  caused  avalanching,  whereas  an  angle  of  33  deg 
seemed  to  be  the  stable  slope  after  avalanching  had  ceased. 

105.  Analysis  of  the  large  wave  tank  data  (Report  1)  showed  smaller 
values  of  the  maximum  slope,  Inferred  to  be  the  angle  of  initial  yield,  with  a 
mean  value  of  28  deg.  The  average  value  for  the  residual  angle  after  shearing 
was  determined  to  be  22  deg;  however,  this  conclusion  is  not  firm  since  some 
profile  steepening  undoubtedly  occurred  between  the  time  of  avalanching  and 
the  next  profile  survey.  Allen  determined  from  experiments  that  the 
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difference  between  4>^y  and  4^  (dilatation  angle)  is  in  the  range  of  10  to 
15  deg  for  sand.  Therefore,  in  the  nuaerical  model,  if  i(>£y  is  set  to  28  deg 
and  4>ra  is  set  to  18  deg,  a  dilatation  angle  of  10  deg  is  used. 

106,  If  the  local  slope  exceeds  4>^y  in  the  model,  material  vlll  be 
redistributed  in  the  neighboring  cells  so  that  the  slope  everywhere  is  leas 
than  or  equal  to  4>ra  .  Computationally  this  is  performed  by  checking  the 
local  slope  at  each  time  step  along  the  grid,  and,  if  4>  >  ^ly  >  avalanch¬ 
ing  is  initiated.  Two  conditions  determine  the  redistribution  of  material; 
mass  should  be  conserved  and  the  residual  angle  after  shearing  is  tra  * 

Figure  6  Is  a  definition  sketch  for  avalanching  where  the  local  slope  between 
cells  1  and  2  exceeds  J’j.y  ■  Depths  before  avalanching  are  denoted  vlth  h  , 
whereas  new  depths  after  sand  redistribution  are  indicated  with  a  prime. 

107.  The  number  of  cells  N  over  which  avalanching  occurs  is  not  known 
a  priori  and  has  to  be  found  iteratively.  Initially,  it  is  assumed  that  the 
avalanching  takes  place  between  cells  1  and  2.  Material  is  moved  from  cell  1 
to  2  in  order  to  obtain  the  slope  4>rl  between  the  cells,  and  the  new  slope 


Figure  6.  Definition  sketch  for  avalanching  calculation 


45 


between  cells  2  end  3  Is  checked  to  see  if  it  exceeds  .  If  this  is  the 

case,  the  calculation  is  repeated,  now  incorporating  cells  1  to  3,  using  the 
mass  conservation  equation  and  assigning  the  slope  fra  to  all  slopes  between 
cells.  This  procedure  is  continued  until  the  slope  between  cells  N  and  N+l 
is  less  than  fra  . 

108.  For  the  general  case  where  sand  is  redistributed  in  H  calcula¬ 
tion  cells  Included  in  avalanching,  che  bsss  conservation  equation  Is  written 

i 

X  Aht  -  0  (39) 

l-i 


where 


Ahj  *  h|  •  hj  (40) 

Here  the  cell  numbering,  denoted  by  1  ,  starts  from  the  cell  where  avalanch¬ 
ing  is  Initiated  and  is  taken  in  the  direction  of  sand  redistribution.  Also, 
depth  changes  have  to  be  taken  with  the  proper  sign.  The  other  criterion 
necessary  to  obtain  the  solution  is  that  after  avalanching,  the  slope  between 
calculation  cells  must  be  fca  .  If  the  length  of  the  calculation  cell  la 
Ax,  the  residual  angle  after  shearing  is  expressed  as 


*ra 


(41) 


Since  fra  Is  constant,  the  difference  in  depth  betveen  neighboring  calls  Ah 
has  to  be  a  constant  given  by 

Ah  -  fra  Ax  (42) 

109.  The  requirement  regarding  the  difference  In  depth  between  cells  as 
given  by  Equation  42  may  be  written  for  N-l  pairs  of  cells  included  in  the 
avalanching  as 
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Ah  —  (hj  +  Ahj)  •  (h,  +  Ah,) 

Ah  -  (h,  +  Ah,)  •  (h,  +  Ah,) 

(43) 

Ah  *>  (hj  +  Aht)  -  (ht+,  +  Ahj.j) 

Ah  -  (hji-i  +  Ah*.,)  -  (h,  +  Ah„) 

This  system  of  equations  may  be  revritten  by  consecutively  eliminating  the 
first  part  of  the  right-hand  side  of  each  aquation  by  adding  the  previous 
equation.  The  result  is 


Ah  -  (hj  +  Ahj)  -  (h,  +  Ah,) 

2Ah  -  (hj  +  Ahj)  -  (h,  ♦  Ah,) 

lAh  -  (hj  +  Ahj)  -  (hj»j  +  Ah,*,) 

(N-l)Ah  -  (hj  +  Ah,)  -  (h*  Ah„) 


(44) 


110.  From  Equation  44  the  unknown  depth  changes  Ah,  may  be  explicitly 
identified  as 


Ah,  *  (hj  +  Ahj)  -  h,  -  Ah 
Ah,  —  (h,  +■  Ah,)  -  h,  -  2Ah 

(45) 

4hj«j  -  (hj  +  Ah,)  -  hj«,  -  lAh 
Ahu  ■*  (h,  +  Ahj)  -  hu  -  (N-l)Ah 

From  Equation  45  it  is  possible  to  compute  all  the  depth  changes  once  the 
depth  change  in  the  first  cell  is  known.  To  determine  the  depth  change  in  the 
first  cell,  the  mass  conservation  equation  (Equation  39)  is  used.  The  sum  of 
ell  changes  in  depth  must  be  zero 
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(46) 


Replacing  the  summation  over  the  depth  change*  from  cell  2  to  N  using 
Equation  45  yield* 


*  *  * 

Ah,  +  l  (h,  +  Ah,)  -  l  h,  -  l  Ah  -  0  (47) 

l-l  l-l  l-l 

111.  From  Equation  47  It  is  possible  to  explicitly  solve  for  the  depth 
change  in  cell  1. 


Ah, 


(N-l)  1 

-  h,  +  - 

N  1  N 


"  1 

l  h,  +  -  (N-l)  Ah 

l-l  2 


(48) 


Depth  changes  In  cells  2  to  N  are  then  given  by 

Ah,  -  h,  +  Ah,  -  h,  -  (i-1)  Ah  (49) 

112.  The  tine  interval  in  the  numerical  model  typically  varies  between 
5  and  20  min  whereas  the  length  step  Is  typically  1  to  5  m.  The  time  interval 
and  length  step  are  not  independent,  since  they  both  govern  the  stability  of 
the  numerical  solution  scheme.  Generally,  a  shorter  length  step  requires  a 
smaller  time  Interval.  An  explicit  criterion  for  the  relation  between  time 
and  length  step  for  numerical  stability  could  not  be  determined  due  to  the 
nonlinearity  of  the  equations  and  the  smoothing  of  the  transport  rata. 

However,  it  vas  empirically  found  that  use  of  the  slope -dependent  term  in  the 
transport  equation  in  the  broken  wave  zone  (Equation  4)  greatly  improves 
stability  of  the  solution.  If  this  term  is  omitted,  numerical  oscillations 
may  originate  from  the  shoreward  side  of  the  bar  as  the  profile  slope  steepens 
under  constant  waves  and  water  level. 


48 


PART  V:  MODEL  TESTS 


113.  Although  the  numerical  model  successfully  reproduced  beach  profile 
change  generated  In  large  wave  tanks.  Including  bar  formation,  bar  movement, 
and  foreshore  erosion,  the  ultimate  test  Is  to  verify  model  performance  for 
field  beaches.  This  chapter  reviews  previous  field  simulations  (Report  1) 
which  concentrated  on  bar  formation  and  movement  and  discusses  recent  field 
simulations  where  the  objective  was  to  test  the  model  for  reproducing  fore¬ 
shore  erosion  as  well  as  accretion.  In  the  field  tests,  documented  storm 
erosion  events  from  both  the  US  east  and  west  coasts  were  simulated.  The 
chapter  begins  with  a  summary  of  sensitivity  analyses  performed  to  investigate 
predictions  of  bar  features  to  changes  In  model  parameters. 

Sensitivity  Analyses  and  Previous  Model  Tests 

114.  The  objective  of  a  sensitivity  analyses  Is  to  establish  the 
Influence  of  model  parameters  on  simulation  results.  Report  1  describes 
simulations  performed  to  quantify  the  effect  of  variations  in  model  parameters 
on  the  prediction.  The  Investigation  focused  on  bar  properties  such  as 
volume,  maximum  height,  and  location  of  the  center  of  mass.  Parameters 
Investigated  were  the  transport  rate  coefficients  K  and  t  ,  breaker  decay 
model  coefficients  k  and  T  ,  and  equilibrium  energy  dissipation  De(]  .  The 
Influence  of  wave  height,  wave  period,  and  runup  height  were  also  explored. 

115.  Crain  size  Is  a  critical  Input  parameter  that  determines  the 
equilibrium  energy  dissipation  and  thus  the  amount  of  material  that  Is  moved 
by  the  waves  in  the  course  of  reaching  an  equilibrium  beach  shape.  Figure  7 
displays  the  time  evolution  of  bar  volume  for  various  grain  sizes  In  a 
hypothetical  eroslonal  case.  The  equilibrium  beach  profile,  defined  as  the 
profile  shape  for  which  no  net  transport  occurs  and  bar  growth  ceases  for  a 
given  wave  condition,  Is  approached  more  rapidly  for  coarser  grain  sizes. 

This  Is  expected  since  less  material  has  to  be  redistributed  by  waves  for  a 
larger  value  of  De(}  .  Grain  size  is  one  of  the  few  model  parameters  that 
significantly  affects  the  movement  of  the  bar  center  of  mass.  Figure  7  Is 
typical  of  the  results  from  the  sensitivity  analysis,  and  similar  curves  are 
presented  in  Report  1  for  all  model  parameters  discussed  here. 
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Figure  7.  Time  evolution  of  ber  volume  for  different 
grain  sizes  for  hypothetical  simulations 


116.  The  transport  rate  coefficient  K  mainly  governs  the  time 
response  of  the  profile  and  has  much  less  Influence  on  bar  volume  and  eroded 
foreshore  sand  volume.  A  large  value  of  K  implies  rapid  profile  response 
with  rapid  attainment  of  equilibrium.  The  coefficient  c  ,  controlling  the 
strength  of  the  slope -dependent  term  In  Equation  4,  does  not  greatly  Influence 
the  time  response,  particularly  at  the  Initial  phase  of  a  simulation,  and 
mainly  contributes  to  reduce  the  depth  of  the  trough  landward  of  the  bar.  Bar 
volume,  however,  Is  affected  by  a  change  In  t  ,  where  a  larger  value  Implies 
a  larger  bar  volume. 

117.  The  empirical  coefficients  in  the  breaker  decay  model  have  a 
marked  influence  on  profile  evolution  with  regard  to  bar  volume.  A  decrease 
In  k  causes  an  increase  in  bar  volume  vhlch  may  be  physically  explained  by 
the  resultant  larger  wave  energy  dissipation  in  the  inner  part  of  the  surf 
zone.  A  smaller  k -value  produces  more  gradual  energy  dissipation,  allowing 
waves  to  retain  more  energy  to  be  dissipated  further  Inshore.  A  large  value 
of  the  stable  wave  height  coefficient  T  Implies  a  flatter  equilibrium  beach 
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profile,  requiring  more  material  to  be  moved  before  equilibrium  Is  attained 
and  thus  more  foreshore  erosion. 

118.  The  effect  of  changing  wave  height  and  wave  period  Is  obvious, 
where  an  Increase  In  either  quantity  Implies  greater  Incident  wave  energy  and 
more  erosion.  Runup  height  has  little  Influence  on  bar  development  since 
about  the  same  amount  of  sand  is  eroded  along  the  profile  Independent  of  runup 
height.  However,  the  shape  of  the  foreshore  can  be  substantially  different 
depending  on  runup  height.  Smaller  runup  height  causes  erosion  to  occur  over 
a  narrower  portion  of  the  foreshore  In  comparison  to  larger  runup;  a  more 
pronounced  scarp  develops  with  a  smaller  runup  height. 

119.  In  Report  1,  the  numerical  model  was  used  to  simulate  profile 
change  measured  at  the  FRF,  located  on  a  barrier  Island  beach  facing  the 
Atlantic  Ocean  at  Duck,  North  Carolina.  High-quality  beach  profile  surveys 
from  the  dune  to  the  depth  of  significant  sand  transport  have  been  made  at 
approximately  two-week  intervals  for  more  than  four  years.  Four  different 
shore-normal  lines  were  monitored  together  with  more  frequent  wave  and  water 
level  measurements.  Cage  measurements  of  wave  height  and  period  are  made  at 
6  hr  intervals,  and  water  level  is  recorded  every  6  min  but  available  at 
hourly  Intervals.  A  portion  of  the  profile  data  set  is  given  in  Howd  and 
Blrkemeler  (1987). 

120.  Five  beach  profile  sets  were  chosen  for  supplying  simulation  data 
for  model  verification.  Each  set  contained  two  profile  surveys  with  a  time 
history  of  waves  and  water  levels  between  surveys.  Erosion  predominated  for 
all  five  events,  including  significant  offshore  movement  of  one  or  two  bars, 
but  shoreline  position  showed  little  or  no  movement.  Shoreline  stability  is 
characteristic  at  the  FRF  beach  and  may  be  caused  In  part  by  an  abundance  of 
coarse  material  on  the  foreshore.  The  data  set  from  Duck  was  thus  unsuitable 
for  testing  predictability  of  the  numerical  model  regarding  shoreline  retreat 
and  subaerial  beach  erosion,  but  the  data  were  considered  appropriate  for 
evaluating  bar  movement. 

121.  Potential  three-dimensionality  of  changes  on  a  natural  beach  must 
be  considered  in  simulating  profile  change.  Even  though  simulation  time 
periods  were  chosen  such  that  neighboring  profile  surveys  most  nearly  dis¬ 
played  similar  development,  some  three-dimensional  effects  were  present.  The 
degree  of  three-dimensionality  was  quantified  by  applying  the  mass  conserve- 
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tion  equation  to  calculate  the  loss  or  gain  of  sand  between  consecutive 
profiles.  Survey  line  188,  located  southeast  of  the  FRF  pier,  was  used  in  the 
simulation  since  It  Is  beyond  major  Influence  of  the  pier,  and  the  bottom 
contours  are  essentially  straight  and  parallel  (Howd  and  Blrkemeier  1987)  . 

122.  Four  of  the  five  events  were  used  for  calibration,  where  the 
calibration  criterion  comprised  minimizing  the  sum  squares  of  the  difference 
between  measured  and  calculated  profile  change.  The  fifth  event  was  used  for 
model  verification  (Report  1) .  The  energy-based  significant  wave  height  (Hmo 
from  gage  620  located  in  18  a  of  water)  was  used  to  calculate  the  transport 
rate  distribution  together  with  the  peak  period  of  the  spectra.  Since  the 
measurement  Intervals  of  waves  and  water  level  were  longer  than  the  model  time 
step,  required  Intermediate  values  were  obtained  by  cubic  spline  interpola¬ 
tion.  A  typical  simulation  time  step  was  20  min,  and  the  length  step  was  5  m. 

123.  Essentially  the  same  parameter  values  obtained  from  the  large  wave 
tank  simulations  were  used,  although  some  modifications  had  to  be  made  for 
field  simulations.  The  most  important  difference  was  that  the  transport 
coefficient  K  had  to  be  decreased  to  obtain  optimal  agreement  for  field 
simulations.  The  optimal  value  of  K  for  field  simulations  was  about  half 
the  value  for  the  large  wave  tank  data.  Two  different  median  grain  diameters 
were  used  to  represent  sediment  characteristics  of  the  FRF  beach  along 
transect  188  (Howd  and  Blrkemeier  1987).  A  median  grain  size  of  2.0  mm  was 
used  on  the  foreshore,  and  0.15  mm  was  used  for  points  seaward  of  the 
foreshore . 

124.  Figure  8a  shows  the  measured  wave  height,  wave  period,  and  water 
level  for  the  simulation  time  period.  Figure  8b  illustrates  one  of  the 
calibration  runs  for  3-6  April  1984  when  a  single  bar  was  present  which  moved 
about  60  m  offshore  during  a  minor  storm  event.  Measured  Initisl  and  final 
profiles  are  given  together  with  the  final  calculated  profile.  Bar  location 
was  rather  well  predicted,  whereas  bar  volume  was  underestimated  and  the 
trough  was  not  sufficiently  pronounced.  Also,  some  shoreline  retreat  was 
calculated  In  spite  of  the  larger  grain  size  on  the  foreshore. 

125.  Model  verification  gave  reasonable  agreement  between  measured  and 
calculated  profiles,  although  the  selected  event  encompassed  profiles  exhibit¬ 
ing  two  bar  features,  which  was  not  the  case  In  the  calibration  events.  As 
with  the  calibration  cases,  direction  of  bar  movement  and  location  were 
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Depth  (m) 


a.  Variation  with  time  of  wave  height, 
period,  and  water  level 


Figure  8.  Calibration  of  numerical  model  against  field  data 
from  Duck,  NC  for  event  840403  •  840406 
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correctly  predicted  whereas  the  size  of  the  bars  was  smaller  than  measured. 
Simulations  performed  subsequent  to  these  early  field  tests  indicate  that  a 
probable  reason  for  the  smaller  calculated  bar  size  Is  Insufficient  varia¬ 
bility  in  the  imposed  wave  Input.  Greater  variability  would  tend  to  shift  the 
break  point  back  and  forth  to  produce  a  wider  bar.  Some  variability  In  wave 
and  water  level  conditions  Is  reproduced  through  use  of  a  time  series  of 
regular  waves,  but  It  is  believed  that  introduction  of  random  fluctuations 
around  these  values  would  give  more  realistic  predictions.  The  absence  of 
randomness  In  the  wave  description  Is  probably  one  of  the  factors  causing  a 
lower  transport  rate  coefficient  in  the  field  simulations.  A  random  variation 
In  wave  height,  for  example,  would  reduce  the  tendency  for  a  steady  trend  in 
profile  evolution,  thereby  requiring  a  larger  value  of  the  transport  rate 
coefficient  K  to  achieve  optimal  agreement. 

126.  Summarizing  experience  from  modeling  beach  profile  evolution  at 
Duck,  prediction  of  the  direction  of  bar  movement  and  bar  location  was 
successful,  although  calculated  bar  size  showed  less  agreement  with  observed 
profile  development.  Profile  response  was  best  reproduced  using  as  much 
information  as  possible  on  wave  height,  wave  period,  and  water  level  (as 
opposed  to  using  average  values  for  the  simulation  period),  lending  credibi¬ 
lity  to  the  predictive  capabilities  of  the  numerical  model.  The  success  of 
the  model  In  predicting  observed  bar  movement  supports  the  hypothesis  that 
breaking  waves  are  the  main  agent  for  cross -shore  sand  movement,  at  least  in 
the  outer  surf  zone.  The  steep  bar  and  foreshore  slopes  appearing  in  the 
large  wave  tank  simulations  were  not  encountered  in  field  simulations  because 
measured  changing  forcing  conditions  were  used.  Reproduction  of  milder  bar 
features  validates  the  concept  of  simulating  profile  change  by  random  waves 
through  a  summation  of  transport  events  produced  by  different  regular  waves 
acting  over  short  time  Intervals. 

Storm  Event  on  the  US  East  Coast 

127.  Several  profile  lines,  located  north  and  south  of  Manasquan  Inlet, 
NJ ,  were  surveyed  by  the  US  Army  Engineer  District,  Philadelphia,  over  a  time 
period  that  encompassed  the  27  to  29  March  northeaster  (Gebert,  Watson,  and 
Hemsley  in  preparation).  The  data  set  includes  profile  surveys  of  the 
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subaerial  beach  taken  1  to  2  days  prior  to  arrival  of  the  storm  and  3-4  days 
after  the  storm.  Pre-storm  subaqueous  profile  data  used  for  simulations  were 
collected  in  early  January  1984  and  post-storm  nearshore  surveys  were  made 
approximately  2  weeks  after  the  event.  A  time  series  of  wave  height  and 
period  was  recorded  by  a  wave  rider  buoy  located  directly  offshore  In  15  m 
(SO  ft)  water  depth;  water  level  at  Kanasquan  Inlet  was  also  available. 

128.  Profile  measurements  were  made  along  8  beach  transects  at 
Manasquan  to  the  north  of  the  Manasquan  Inlet  jetties  and  along  9  transects  at 
Point  Pleasant  Beach  to  the  south  of  the  Jetties  (Figure  9).  Point  Pleasant 
Beach  Is  relatively  wide  and  sandy  with  no  coastal  structures  in  the  surf  zone 
south  of  the  678-m  long  south  jetty  at  Manasquan  Inlet,  whereas  Manasquan  has 
a  more  narrow  sandy  beach  with  short  groins  located  at  approximately  200-  to 
300-m  intervals. 

129.  Regular  profile  surveys  were  carried  out  at  Kanasquan  and  Point 
Pleasant  Beach  from  1982  to  1984  under  the  auspices  of  the  Corps  of  Engineers 
Monitoring  Completed  Coastal  Projects  (MCCP)  Program.  The  main  purpose  of  the 
project  was  to  monitor  the  structural  stability  of  the  Jetties  after  rehabili¬ 
tation  with  16-ton  dolosse  armor  units.  Two  secondary  concerns  dealt  with 
maintenance  of  inlet  cross-sectional  area  and  potential  change  on  adjacent 
beaches  caused  by  the  Jetties  (Gebert,  Watson,  and  Hemsley  in  preparation). 
Figure  9  gives  the  locations  of  the  survey  lines  and  the  numbering  convention. 
The  lines  are  numbered  consecutively  from  1  starting  at  the  Jetties  and 
increasing  north  and  south,  where  the  symbols  M  and  PP  are  used  to  denote 
Manasquan  and  Point  Pleasant  Beach,  respectively.  A  tide  gage,  operated  by 
the  National  Oceanographic  Service  (NOS),  is  located  on  the  bay  side  of 
Manasquan  Inlet.  At  the  time,  a  wave  rider  buoy  was  operated  by  MCCP  and 
located  in  15  m  (50  ft)  of  water  seaward  of  Manasquan  (Figure  9).  Wave 
conditions  generally  were  measured  every  6  hr  for  a  20-mln  time  period,  and 
the  sampling  interval  was  automatically  decreased  to  1  or  2  hr  when  the  wave 
height  exceeded  3  m. 

130.  Toward  the  end  of  March  1984,  a  northeaster  passed  the  area 
causing  severe  erosion  of  the  beaches.  The  time  history  of  measured  signifi¬ 
cant  wave  height  and  peak  spectral  wave  period  for  the  storm  is  displayed  in 
Figure  10  and  water  level  is  shown  in  Figure  11.  Water  surface  elevation  is 
given  with  respect  to  mean  sea  level  (MSL)  and  comprises  both  storm  surge  and 
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Figure  9.  Location  map  shoving  profile  survey  lines  and  the  tide 
gage  and  wave  rider  buoy  positions  at  the  Manasquan  and 
Point  Pleasant  Beach  study  area 


tidal  variations.  Wave  height  began  to  increase  In  the  morning  on  28  Harch 
and  peaked  about  24  hr  later  on  29  Harch  with  a  maximum  height  of  6.6  m.  By 
noon  on  30  March  the  storm  had  passed  and  more  typical  wave  heights  were 
recorded.  Have  period  Increased  sharply  vlth  arrival  of  the  storm,  similar  to 
the  behavior  of  wave  height,  but  long  periods  persisted  after  the  storm 
creating  mlld-steepness  waves  and  strongly  accretionary  conditions. 


Wave  Height  and  Period 


Figure  10.  Uave  height  and  period  measured 
northeast  of  Kanasquan  Inlet 

Water  Level  (m) 


Figure  11.  Water  level  measured  on  the  landward 
side  of  Manasquan  Inlet 
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131.  The  water  level  record  displayed  a  pronounced  peak  for  the  middle 
of  the  storm,  approximately  coinciding  with  the  maximum  wave  height.  The 
tidal  variation  at  Manasquan  Inlet  Is  semidiurnal,  with  two  high  and  two  lov 
waters  in  a  tidal  day,  and  the  mean  tidal  range  is  1.2  m.  Maximum  surface 
elevation  measured  during  the  storm  event  was  1.90  m  at  about  0600  on  29  March. 
Profile  measurements 

132.  Pre-storm  beach  profile  surveys  were  taken  on  26  and  27  March  for 
Manasquan  and  Point  Pleasant  beach,  respectively.  However,  when  the  storm 
arrived,  the  scheduled  fathometer  surveys  had  to  be  halted  due  to  high  waves, 
and  only  wading  depth  profiles  were  obtained  for  these  dates.  Subaqueous 
profile  measurements  obtained  from  the  previous  profile  survey  (6  January  1984 
at  Manasquan  and  18  January  1984  at  Point  Pleasant  Beach)  were  used  to 
construct  complete  pre-storm  profiles.  Post-storm  profiles  were  recorded  on 

3  April  (onshore)  and  11  April  (offshore)  at  Manasquan  and  2  April  (onshore) 
and  12  April  (offshore)  at  Point  Pleasant  Beach.  Profile  survey  elevations 
were  made  with  respect  to  NGVD  (National  Geodetic  Vertical  Datum)  ,  and  were 
converted  to  mean  sea  level  In  this  report  by  the  relationship  MSL(m)  - 
NCVD(m)  +  0.24  m. 

133.  Profile  survey  data  were  analyzed  to  yield  representative  profiles 
to  be  used  In  numerical  simulations.  Also,  longshore  variability  In  eroded 
volume  was  Investigated  to  evaluate  the  applicability  of  a  two-dimensional 
approach.  Since  the  shoreline  was  curved  close  to  the  jetties,  especially  at 
Point  Pleasant  Beach,  profiles  at  the  respective  beaches  were  shifted  across- 
shore  and  aligned  with  reference  to  MSL.  A  single  representative  profile  for 
each  beach  was  determined  by  averaging  profile  elevations  at  fixed  Intervals. 

134.  Figure  12  shows  plots  of  the  eight  pre-storm  profiles  measured  at 
Point  Pleasant  Beach.  The  shapes  of  the  profiles  are  surprisingly  similar, 
Indicating  relatively  good  two -dimensionality.  However,  the  two  profiles 
located  closest  to  the  jetty  (PP-1  and  PF-2)  show  depths  that  are  signifi¬ 
cantly  less  than  those  associated  with  the  other  profiles.  This  shallower 
subaqueous  portion  of  the  profile  Is  probably  a  result  of  Impoundment  of  sand 
transported  alongshore.  The  two  profiles  closest  to  the  jetty  were  therefore 
excluded  when  determining  the  generic  profile  for  simulation.  These  two 
profiles  also  exhibited  markedly  less  erosion  during  the  storm  as  compared 
with  the  other  profiles. 
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Elevation  (m) 


Figure  12.  Profiles  measured  at  Point  Pleasant  Beach  during 
January  1984  aligned  vlth  respect  to  MSL 

135.  In  a  manner  similar  to  the  analysis  for  Point  Pleasant  Beach,  two 
profiles  were  removed  from  the  analysis  at  Manasquan;  seven  profiles  were  used 
for  determining  a  generic  profile  for  the  simulation.  The  transect  closest  to 
the  jetty  (M-l)  was  not  used  along  with  the  profile  furthest  to  the  north 
(H-9).  Their  exclusion  was  motivated  by  the  significantly  different  response 
of  the  two  profiles  to  the  storm,  for  which  substantial  accumulation  took 
place  at  their  seaward  ends.  This  is  attributed  to  longshore  effects  and  the 
Influence  of  structures. 

136.  A  generic  pre-storm  profile  was  created  using  the  average  profile 
from  the  immediate  pre- storm  subaerial  surveys  matched  with  an  average 
subaqueous  profile  developed  from  January  1984  survey  data.  In  the  same 
manner,  a  post-storm  profile  was  developed  based  on  the  surveys  of  2  and  3 
April  and  11  and  12  April  1984. 

137.  Representative  pre-  and  post-storm  profiles  determined  for 
Manasquan  and  Point  Pleasant  Beach  (Figures  13  and  14)  differ  in  two  charac- 
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teristics.  First,  the  dune  Is  higher  end  narrower  at  Manasquan  as  compared  to 
Point  Pleasant  Beach.  Second,  a  wide  subaqueous  bar  exists  In  the  profile  at 
Manasquan,  extending  approximately  300  m  seaward  of  the  baseline,  out  to  the 
end  of  the  groins.  It  appears  as  If  the  groin  field  Is  effective  In  trapping 

-  sand  on  the  nearshore  portion  of  the  profile. 

Etath  tialn.jiiz.t 

138.  Grain  size  is  a  principal  parameter  in  the  numerical  model  and  can 
be  difficult  to  quantify  for  use  In  SBEACH  since  a  wide  range  of  grain  sizes 

-Is  normally  present  on  a  natural  beach.  In  the  study  area,  investigations 
have  been  made  on  the  temporal  and  spatial  distribution  of  grain  size  (Corps 
of  Engineers  1954;  Gorman,  Stauble,  and  Reed  in  press).  However,  because  of 
inconvenient  sample  location  and  lack  of  measurements  in  the  dune  region, 
available  grain  size  information  was  supplemented  with  sampling  from  the 
subaerial  portion  of  the  beach. 

139.  The  additional  sediment  sampling,  performed  by  Mr.  Jeffrey  Gebert, 
US  Army  Engineer  District,  Philadelphia,  in  April  1989,  showed  a  uniform 
median  grain  size  of  around  0.5  mm  In  the  dune  region.  Median  grain  size 
varied  more  on  other  portions  of  the  subaerial  profile  with  coarser  material 
typically  present  at  the  berm  crest.  Sand  samples  were  taken  at  four  cross- 
shore  locations  on  the  dry  portion  of  the  beach  and  at  four  alongshore  sites. 
Table  1  gives  a  summary  of  median  grain  size  data  available  for  the  Manasquan 
and  Point  Pleasant  Beach  study  site. 

140.  Two  sampling  sites  south  of  the  Jetties  and  one  to  the  north  were 
used  in  defining  median  grain  size  In  the  dune  region,  together  with  data  from 
previous  investigations .  Samples  collected  adjacent  to  the  Jetties  were 
excluded  from  the  analysis  by  reasoning  similar  to  chat  used  In  determining  a 
representative  profile.  An  average  median  grain  size  for  the  dune  region  was 
determined  to  be  0.49  mm  for  Manasquan  and  0.54  mm  for  Point  Pleasant  Beach. 

On  the  foreshore,  coarser  material  was  found  near  the  berm  crest  with  median 
grain  sizes  around  0.3  to  1.5  mm.  For  the  subaqueous  portion  of  the  profile, 
median  grain  diameter  from  previously  collected  sand  samples  showed  an  overall 
trend  of  decreasing  grain  size  with  distance  offshore.  The  average  median 
grain  size  for  all  subaqueous  samples  showed  Djq  -  0.19  mm  for  Manasquan  and 
D50  -0.22  mm  for  Point  Pleasant  Beach. 


60 


Elevation  (m) 


Figure  13.  Generic  pre-  and  post-storn  profile  for  Manasquan 
Elevation  (m) 


Pre-storm  - Post-storm 

Figure  1G.  Generic  pre-  and  post-storn  profile  for  Point  Pleasant  Beach 
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Table  1 

Median  Grain  Size  (urn)  for  Manasauan  and  Point  Pleasant  Beach 


Landward  Limit 


of  Dune 

Dune  Crest 

Berm  Crest 

Mid-Tide 

Nearshore3 

Manasauan 

M-l2 

0.38 

0.43 

0.20 

M-2* 

.... 

0.47 

0.56 

0.86 

.... 

M-92 

.... 

.... 

0.44 

0.38 

0.16 

M-93 

.... 

0.51 

0.60 

1.31 

0.21 

Point  Pleasant 

PP-12 

Beach 

0.28 

0.34 

0.21 

PP-1J 

.... 

0.53 

.... 

0.36 

0.29 

PP-1* 

0.54 

0.55 

1.51 

0.51 

.... 

PP-3* 

0.47 

0.53 

1.13 

0.82 

.... 

PP-6* 

0.48 

.... 

1.05 

0.68 

.... 

PP-82 

.... 

.... 

0.53 

0.64 

0.19 

PP-83 

.... 

0.66 

0.41 

0.47 

0.18 

1  average  DJ0  on  subaqueous  profile;  2  June  1953;  *  May  1987;  *  April  1989 

141.  An  Indirect  method  of  obtaining  a  representative  grain  size  for 
the  subaqueous  portion  of  the  profile  Is  to  approximate  the  beach  profile  with 
the  equation  proposed  by  Dean  (1977)  and  solve  for  the  shape  parameter.  From 
the  shape  parameter  It  Is  possible  to  determine  the  corresponding  grain  size 
based  on  curves  presented  by  Moore  (1982)  and  Dean  (1987).  The  equation 
h  -  Ax^/^  proposed  by  Dean  (1977)  for  describing  the  equilibrium  beach 
profile  vas  least-square  fitted  to  the  pre-storm  profile  from  the  January 
surveys  to  Indirectly  determine  a  representative  grain  size  for  the  subaqueous 
portion  of  the  profile.  In  this  equation,  the  quantity  A  is  an  empirical 
function  of  grain  size.  The  effective  median  grain  diameter  calculated  for 
profiles  M-2  to  M-8  was  0.47  mm,  and  for  profiles  PP-3  to  PP-8  the  effective 
grain  size  was  0.51  mm.  In  the  least-square  fit,  both  the  shape  parameter  and 
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the  horizontal  location  of  the  still-water  shoreline  were  estimated  to  minimize 
the  sum  of  squares  between  measured  profiles  and  the  profile  equation. 

142.  The  profile  equation  could  in  general  explain  a  large  portion  of 
the  variation  in  the  measurements  but  visual  inspection  showed  some  deficien¬ 
cies  in  the  description.  To  illustrate  this  point,  Figure  15  shows  the 
representative  pre-storm  profile  at  Point  Pleasant  Beach  and  the  least-square 
fitted  profile.  The  corresponding  median  grain  size  for  the  least-square 
fitted  profile  Is  0.5  mm.  The  fitted  profile  gives  an  overall  reasonable 
description  but  deviates  significantly  compared  to  the  representative  profile. 
Directly  seaward  of  the  shoreline,  the  representative  profile  Is  steeper  than 
the  fitted  profile  whereas  the  opposite  occurs  In  the  offshore  region.  The 
actual  profile  shape  is  probably  sensitive  to  varying  median  grain  size, 
resulting  In  a  shape  which  locally  deviates  from  a  fit  involving  a  2/5  power 
for  a  single  grain  size. 

Alongshore  variability  in  storm  erosion 

143.  The  numerical  model  simulates  cross-shore  sand  transport  only; 
longshore  effects  are  not  represented.  In  order  to  evaluate  this  approxi¬ 
mation  for  the  target  storm  event,  the  amount  of  erosion  was  determined  for 
each  profile  survey.  Equal  amounts  of  erosion  would  suggest  that  the  forcing 
conditions  were  similar  during  the  storm  event  and  differentials  In  longshore 
transport  were  small.  It  is  important  to  note  that  significant  recovery 
occurred  before  the  post-storm  profile  was  surveyed.  Figure  16  Illustrates 
the  recorded  volume  difference  above  the  1-m  contour  between  the  pre-  and 
post- storm  profile.  Profile  location  is  referenced  to  the  profile  closest  to 
the  jetties,  going  north  and  south  at  Manasquan  and  Point  Pleasant  Beach, 
respectively.  Table  2  lists  measured  volumes  of  erosion  above  the  1-m  and  0-m 
(MSL  shoreline)  contours. 

144.  Significantly  less  erosion  is  Indicated  for  the  two  profile  lines 
closest  to  the  Jetties  at  Point  Pleasant  Beach.  South  of  profile  line  PP-2 
the  eroded  volume  Is  fairly  constant  (50  to  65  mV*)  except  for  line  PP-6. 
Examination  of  the  post-storm  profile  survey  at  line  PP-6  shows  considerable 
deposition  on  top  of  the  dune.  As  previously  noted,  the  eroded  volume 
calculation  Interval  includes  a  period  of  recovery  or  accretion  that  contami¬ 
nates  the  data  set  with  respect  to  storm- induced  erosion  modeling.  The 
average  volume  eroded  above  the  1-m  contour  for  all  profile  lines  at  Point 
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Elevation  (m) 


Pre-atorm  - Least-square 

Figure  15.  Generic  pre-storm  profile  at  Point  Pleasant  Beach 
and  least-square  fitted  equilibrium  profile 

Eroded  Volume  (m»»3/m) 


— —  Point  Pleasant  Beach  Manaaquan  i 

Figure  16.  Volume  difference  above  the  1-m  contour  between  pre- 
and  post-storm  profiles  at  Manasquan  and  Point  Pleasant  Beach 
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Table  2 


Eroded  Profile  Volumes  at  Manasquan  and  Point  Pleasant  Beach 
for  the  March  1984  Storm  Event  (nVi») 


Groun  Averaee 

Average  of  Selected 
Profiles 

Minimum 

Maximum 

Manasquan 

0- Contour 

36.3 

39.3 

14.7 

62.6 

1-m  Contour 

32.4 

34.7 

13.4 

56.2 

Point  Pleasant  Beach 

0- Contour 

48.9 

55.8 

9.6 

66.0 

1-m  Contour 

45.8 

53.8 

8.5 

56.2 

Pleasant  Beach  Is  45.8  m3/m,  and  for  profiles  PP-3  to  PP-8  the  average  eroded 
volume  is  53.8  m3/m.  To  gage  the  intensity  of  storm  activity  and  resultant 
magnitude  of  beach  profile  change,  data  from  this  event  were  compared  with 
similar  parameters  from  13  other  East  Coast  storms  (Blrkemeier,  Savage,  and 
Leffler  1988).  Comparison  of  peak  wave  and  water  level  conditions,  and 
associated  changes  in  beach  profile  volume,  with  data  from  the  27  to  29  March 
1984  event  at  Manasquan  and  Point  Pleasant  Beach  indicate  that  dynamic  charac¬ 
teristics  associated  with  this  storm  and  resultant  volume  changes  were  more 
extreme  than  any  of  the  13  previously  recorded  events.  Analysis  of  extreme 
wave  statistics  for  this  area  suggests  that  this  major  storm  event  has  a 
frequency  of  occurrence  of  about  30  years  (Jensen  1983,  Corson  and  Tracy  1985). 

145.  The  profiles  at  Manasquan  experienced  less  erosion  than  those  at 
Point  Pleasant,  with  an  average  eroded  volume  for  all  profile  lines  of  32.4 
m3/m.  Survey  line  M-l  closest  to  the  jetty  experienced  somewhat  less  erosion 
than  all  but  one  of  the  other  profiles.  The  average  volume  eroded  above  the 
1-m  contour  for  profiles  M-2  to  M-8  was  34.7  m3/m. 

146.  Variability  in  volume  of  erosion  is  greater  for  Manasquan  than  for 
Point  Pleasant  if  the  profile  lines  closest  to  the  jetty  are  neglected  (Figure 
16).  This  is  probably  due  to  localized  processes  associated  with  the  groin 
field,  giving  rise  to  longshore  variability  in  sand  transport.  The  lower 
average  eroded  volume  for  Manasquan  is  also  partially  attributed  to  the 
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groins,  which  evidently  caused  a  large  subaqueous  bar  to  develop.  This 
feature  increases  energy  dissipation  on  the  seaward  portion  of  the  profile  and 
reduces  wave  energy  reaching  the  dune  face.  Expressed  differently,  the  bar  is 
a  sand  reservoir  which  must  first  be  depleted  by  erosive  waves  before  they  can 
fully  attack  the  subaerial  beach. 

147.  In  summary,  it  appears  that  longshore  variability  in  eroded  volume 
was  relatively  small  for  Point  Pleasant  Beach  because  structures  were  not 
present  to  create  differentials  in  longshore  transport.  At  Manasquan,  larger 

-volume  differences  are  noted  due  to  the  presence  of  structures  and  the 
location  of  profile  lines  with  respect  to  the  groins.  Applying  the  numerical 
model  to  a  representative  generic  profile  should  be  acceptable  both  at 
Manasquan  and  Point  Pleasant  Beach  in  terms  of  longshore  variability  in  eroded 
volume.  However,  if  only  one  profile  survey  were  available  for  modeling, 
significant  variability  or  deviation  in  calculated  and  measured  results  could 
be  expected. 

Numerical  simulations 

148.  Average  representative  profiles  were  used  in  the  numerical  simula¬ 
tions.  Since  median  grain  diameter  did  not  differ  significantly  between  the 
two  sites,  the  same  value  was  specified.  A  schematized  description  of  cross¬ 
shore  variation  in  grain  size  was  specified  based  on  sampling  results,  with  a 
diameter  of  0.5  mm  shoreward  of  the  MSL  shoreline  and  0.3  mm  seaward  of  the 
MSL  shoreline.  Model  parameters  were  initially  set  at  the  same  values  as  for 
the  large  wave  tank  studies  but  some  modifications  were  made  in  the  calibra¬ 
tion  process. 

149.  Simulation  periods  for  the  two  beaches  corresponded  with  different 
survey  times  at  the  respective  sites.  The  exact  time  of  day  of  the  surveys 
was  not  known  and  was  assumed  to  be  1200.  For  Manasquan  the  simulation 
Interval  was  26  March  to  3  April,  and  for  Point  Pleasant  Beach  it  was  27  March 
to  2  April.  Wave  and  water  level  conditions  for  the  calculation  periods  are 
given  in  Figures  10  and  11,  for  which  cubic  spline  interpolation  was  used  to 
obtain  intermediate  values  between  measurement  points.  The  model  time  step 
was  5  min,  and  the  length  step  was  2  m. 

150.  Simulation  results  are  shown  in  Figure  17  for  Manasquan  and  in 
Figure  18  for  Point  Pleasant  Beach.  The  initial  and  final  representative 
profiles  are  plotted  with  the  calculated  post-storm  and  final  profile.  For 
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Figure  17.  Measured  and  calculated  profile  change  at  Manasquan 


Elevation  (m) 


Figure  18. 


Measured  and  calculated  profile  change  at 
Point  Pleasant  Beach 
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the  calculation,  "post- storm"  was  defined  as  the  time  at  which  the  model 
switched  from  calculating  erosion  to  accretion  by  Equation  1,  and  "final"  was 
defined  as  the  time  at  which  the  post-storm  survey  was  taken.  Since  post¬ 
storm  profile  measurements  were  made  after  several  days  of  recovery  waves,  the 
final  measured  profile  Includes  a  period  of  potentially  significant  beach 
recovery.  In  particular,  at  Point  Pleasant  the  recovery  was  very  strong  with 
a  wide  gentle  ridge  evolving  with  a  crest  2  m  above  the  MSL  shoreline.  Beach 
recovery  at  Manasquan  was  considerable  but  not  as  pronounced  as  at  Point 
Pleasant  Beach,  probably  due  in  part  to  a  steeper  dune  face  and  a  narrower 
foreshore,  conditions  which  tend  to  retard  the  recovery  process. 

151.  Manasquan.  Simulated  profile  retreat  at  the  top  of  the  dune  is 
close  to  the  measured  retreat  for  Manasquan,  and  the  dune  slope  is  reproduced 
well.  Also,  the  post-storm  eroded  volume  is  considered  to  be  in  good  agree¬ 
ment  with  what  probably  occurred  as  inferred  from  the  profile  taken  on  3  April 
by  extrapolating  the  slope  of  the  dune  face  to  MSI  and  neglecting  the  accre¬ 
tionary  berm  feature.  Model  simulation  of  recovery  is  limited  to  qualitative 
and  not  quantitative  agreement  with  the  measurements,  showing  berm  build-up 
but  with  a  smaller  volume.  Another  aspect  of  profile  development  which 
differs  between  calculated  and  measured  is  the  evolution  of  the  bar  located 
near  the  seaward  terminus  of  the  groins.  The  numerical  model  predicts  that 
the  upper  part  of  the  terrace  is  transported  further  offshore  whereas  measure¬ 
ments  indicated  that  this  feature  was  stable  during  the  event.  This  short¬ 
coming  of  the  model  may  be  partially  attributed  to  neglect  of  longshore 
processes . 

152.  Point  Pleasant  Beach.  The  model  simulation  for  Point  Pleasant 
Beach  shows  much  the  same  characteristics  as  for  Manasquan  although  predicted 
dune  retreat  was  smaller  as  compared  to  measured  retreat.  Post-storm  eroded 
volume  is  again  considered  to  be  in  good  agreement  with  the  actual  amount 
through  examination  of  the  2  April  profile  and  disregarding  the  large  ridge 
which  was  assumed  to  have  been  formed  during  the  recovery  phase.  Simulated 
recovery  is  stronger  at  Point  Pleasant  compared  to  Manasquan  which  is  qualita¬ 
tively  in  accordance  with  observations,  although  the  predicted  berm  volume  is 
smaller.  The  numerical  model  also  indicates  development  of  bar-like  features 
in  the  seaward  portion  of  the  profile  which  were  not  present  in  the  measured 
profile . 
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Discussion 

153.  In  the  simulations  for  both  beaches,  the  transport  rate  coeffi¬ 
cient  K  was  changed  to  obtain  optimal  agreement  between  model  results  and 
the  measured  profile.  The  coefficient  value  determined  for  both  sites  was 

K  -  2.0  10"6  m^/N  ,  compared  with  the  value  from  the  large  wave  tank  simula¬ 
tions  of  K  -  1.6  10‘®  n/*/N  and  from  the  Duck,  North  Carolina,  simulations  of 
K  -  0.7  10' ®  a*/H  .  The  Manasquan  and  Point  Pleasant  Beach  profile  change 
simulations  used  an  optimal  K-value  similar  to  the  large  wave  tank  simulation, 
whereas  the  coefficient  for  the  Duck  data  was  considerably  lower,  probably  due 
to  the  unusual  stability  of  the  shoreline  at  Duck  that  hindered  erosion. 

154.  To  better  account  for  naturally  variable  wave  conditions,  wave 
heights  used  in  model  simulations  were  randomly  varied  a  certain  percentage 
around  the  values  interpolated  from  the  measured  wave  heights.  This  technique 
simulates  variability  in  the  location  of  the  break  point  and  development  of 
profile  features.  The  method  for  randomly  varying  wave  height  consisted  of 
choosing  a  height  from  a  rectangular  distribution  with  its  center  as  the 
interpolated  height  and  its  width  as  a  predefined  percentage  of  the  interpo¬ 
lated  height.  Different  distribution  widths  were  tested,  and  10  to  20  percent 
proved  to  give  a  good  qualitative  description  of  profile  evolution.  For  the 
calibration  runs  displayed  in  Figures  17  and  18,  a  20  percent  variation  was 
employed.  Although  this  technique  was  not  performed  in  simulations  for  the 
case  at  Duck,  a  larger  K-value  and  smoother,  more  realistic  profile  features 
would  be  expected. 

155.  In  order  to  obtain  a  qualitatively  improved  description  of  the 
accretion  process,  adjustments  to  the  transport  rate  distribution  were  made. 
For  the  large  wave  tank  data  simulations,  the  transport  rate  distribution  was 
calculated  the  same  way,  independent  of  whether  onshore  or  offshore  transport 
occurred.  However,  for  a  natural  beach  with  a  variable  wave  climate,  unusual 
conditions  arise  which  must  be  dealt  with  in  a  special  way.  For  example,  some 
accretionary- type  waves  are  small  enough  such  that  breaking  will  not  occur 
until  they  arrive  almost  at  the  shoreline.  The  transport  rate  distribution  in 
this  case  is  very  narrow  and  peaked,  producing  large  accretion  in  the  cells 
around  the  shoreline.  This  local  accumulation  steepens  the  lower  foreshore 
slope,  making  more  waves  break  at  the  shoreline.  To  take  care  of  this  condi¬ 
tion,  if  at  any  time  step  the  surf  zone  was  calculated  to  be  less  than  four 
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cells  wide,  che  transport  rate  was  specified  as  linearly  decreasing  from  the 
plunge  point  to  the  runup  limit.  The  transport  rate  at  the  plunge  point  was 
determined  from  the  wave  energy  dissipation  according  to  Equation  4,  as  usual. 

156.  If  onshore  transport  is  high,  berm  build-up  Is  rapid  and  the 
foreshore  beach  slope  may  become  unrealistically  steep,  United  only  by  the 
angle  of  Initial  yield.  Such  slopes  predicted  by  the  nunerlcal  model  are  In 
contradiction  with  observations  from  the  large  wave  tank  data  (Takeda  and 
Sunamura  1983,  Larson  and  Kraus  1989a)  and  much  steeper  than  what  la  found  on 
a  natural  beach  (Bascom  1951,  Ulegel  1964,  Thomas  and  Baba  1986,  Hughes  and 
Cowell  1987).  To  limit  steepening  of  the  foreshore,  the  cross-shore  sand 
transport  rate  was  reduced  with  respect  to  the  local  slope  when  transport  was 
directed  onshore.  This  Introduced  a  second  dependence  on  beach  slope  since, 
Equation  4  already  contains  a  term  dependent  on  dh/dx  ,  but  was  necessary  to 
achieve  a  smoother,  more  realistic  approach  to  equilibrium. 

157.  Eroded  volumes  for  the  predicted  post-storm  profile  and  profile 
after  recovery  were  calculated  above  the  0-  and  1-m  contour  and  are  displayed 
In  Table  3.  These  values  can  be  compared  with  eroded  volumes  based  on 
measured  profile  evolution  listed  in  Table  2.  As  previously  discussed,  eroded 
volume  determined  from  the  measured  profiles  Includes  a  period  of  recovery  as 
the  Immediate  post-storm  profile  was  not  surveyed.  Since  model  prediction  of 
recovery  was  only  In  qualitative  agreement  with  measurements,  calculated  eroded 

Table  3 

Calculated  Eroded  Profile  Volumes  at  Manasauan 
and  Point  Pleasant  Beach  U aVm) 


Maximum  Erosion 

With  Recovery 

Manas quan 

0- Contour 

86.4 

73.9 

1-m  Contour 

54.6 

54.2 

Point  Pleasant  Beach 

0- Contour 

51.5 

44.2 

1-m  Contour 

37.2 

37.2 
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volumes  at  the  end  of  the  simulation  period  were  too  large  In  comparison  with 
measured  changes.  Calculated  eroded  volume  at  Manasquan  agrees  well  with 
measurements,  although  the  shape  of  the  profile  differs,  with  the  measured 
profile  exhibiting  less  foreshore  erosion  and  greater  dune  retreat.  Post¬ 
storm  erosion  Is  probably  well  predicted  at  Point  Pleasant  Beach  but  berm 
bulld-up  is  underestimated.  Little  or  no  recovery  Is  produced  by  the  numeri¬ 
cal  model  above  the  1-m  contour  which  differs  from  what  occurred  in  the  field. 

158.  In  conclusion,  the  numerical  model  was  successful  in  reproducing 
the  erosional  portion  of  the  studied  storm  event  but  predicted  profile 

_ evolution  during  the  recovery  stage  was  only  In  qualitative  agreement  with 
measurements.  Parameter  values  obtained  for  model  calibration  were  similar  to 
those  found  for  simulations  of  the  large  wave  tank  data,  particularly  the 
transport  rate  coefficient  K  .  A  two-dimensional  approach  to  estimate 
profile  evolution  during  a  storm  event  was  Justified  based  on  the  small 
measured  alongshore  variability  in  eroded  volume.  However,  if  nearshore 
structures  such  as  groins  exist  in  a  study  area,  care  should  be  taken  in 
choosing  a  representative  profile,  and,  if  possible,  averages  of  several 
profiles  should  be  used. 

Sti»rm  Evgqt_pn_iht-!Jg  West  Cpaat 

159.  A  mission-support  project  performed  by  CERC  for  the  Coastal  of 
California  Storm  and  Tidal  Waves  Study  conducted  by  US  Army  Engineer  District, 
Los  Angeles,  provided  an  opportunity  to  apply  SBEACH  to  storm-erosion  events 
on  the  southern  coast  of  California.  This  application  involved  considerably 
different  environmental  conditions  than  East  coast  field  tests  of  the  model. 
West  coast  conditions  of  note  include  small  storm  surge,  typically  less  than 
0.3  m,  longer-period  waves  than  on  the  East  coast,  and  broad,  low-relief  bars 
if  bars  are  present. 

160.  Documentation  of  storm- induced  beach  erosion  suitable  for  model 
testing  and  calibration  is  relatively  scarce  for  the  West  coast.  As  of  late 
1988,  after  a  search  of  the  literature  and  inquiry  to  leading  coastal  investi¬ 
gators,  only  the  data  collected  by  Aubrey  (1978)  was  considered  to  be  adequate. 
The  measurements  of  interest  were  made  at  Torrey  Pines  Beach,  California, 
located  approximately  2  miles  north  of  Scrlpps  Institution  of  Oceanography 
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(Aubrey,  Inman,  and  Nordstrom  1976;  Aubrey  1978;  Gable  1979),  between  December 
1975  and  January  1977. 

161.  Figure  19  shows  the  locations  of  three  range  lines  along  Torrey 
Pines  Beach  referred  to  as  North  Range,  Indian  Canyon  Range,  and  South  Range. 
Onshore  and  offshore  portions  of  the  profiles  were  measured  periodically 
between  1972  and  1976,  together  with  the  associated  wave  climate  and  water 
level.  All  elevations  were  measured  relative  to  NGVD;  mean  sea  level  la 

4.3  cm  above  NGVD  In  this  region.  A  time  series  of  waves  and  water  level  for 
the  periods  29  December  1975  to  9  January  1976  and  21  December  1976  to  4 
January  1977  were  available  for  defining  short-term  beach  profile  change  due 
to  storm  conditions.  Vave  measurements  were  collected  along  Indian  Canyon 
Range  with  an  array  of  pressure  sensors  located  at  a  depth  of  approximately 
9.1m  (Figure  19).  Examination  of  the  Torrey  Pines  data  shows  these  events  to 
be  less  significant  than  most  storms  producing  notable  beach  change  In 
southern  California.  However,  this  Information  provides  the  only  compre¬ 
hensive  data  base  on  storm  conditions  and  associated  beach  profile  response 
for  the  region  that  could  be  found  during  this  study. 

162.  Beach  profile  data  were  collected  three  times  during  the  first 
storm  event  (12/29/75,  12/31/75,  and  1/9/76)  and  twice  during  the  second  event 
(12/21/76  and  1/4/77).  All  surveys  were  carried  out  near  the  time  of  predict¬ 
ed  low  water.  The  offshore  bathymetry  indicated  very  little  longshore 
variation,  Implying  only  minor  differences  In  wave  properties  at  the  three 
range  lines.  Therefore,  an  average  profile  shape  was  constructed  for  each 
date  by  matching  the  zero  crossings  (referenced  to  NGVD)  for  each  of  the  three 
profiles  and  averaging  the  elevations  at  known  positions  relative  to  y  -  0. 
Use  of  an  average  profile  should  reduce  changes  associated  with  longshore 
variability  caused  by  local  differences  in  nearshore  bottom  topography, 
sediment  size,  and  waves  (Savage  and  Blrkemeier  1987,  Blrkemeler  et  al.  1987). 
These  composite  data  were  used  for  simulating  beach  profile  response  to  both 
storm  events.  The  recorded  time  series  of  waves  and  water  level  between 
surveys  was  Input  to  the  beach  profile  change  model  SBEACH,  and  data  quality 
checks  were  made  on  all  process  and  response  measurements  before  model  runs 
proceeded. 

163.  Model  parameter  values  resulting  from  calibration  with  the  large- 
scale  laboratory  data  sets  (Report  1)  were  Initially  used  to  simulate  field 
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Figure  19.  Location  map  of  Torrey  Pines  Beach  showing  position  of  range 
lines  and  pressure  sensor  array  (after  Aubrey  1978) 

beach  profile  response.  It  became  apparent,  however,  that  values  of  the  two 
principal  empirical  coefficients  would  have  to  be  modified  to  achieve  accept¬ 
able  agreement  betveen  measured  and  calculated  profiles.  The  optimal  sand 
transport  rate  coefficient  K  In  Equation  4  obtained  from  these  two  West 
coast  storm  events  was  slightly  smaller  than  that  obtained  from  the  laboratory 
experiments  (1.4x10'®  m^/N  versus  1.6x10'®  m^/N) ,  and  the  transport  rate 
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coefficient  t  for  the  slope-dependent  tern  of  the  transport  rate  equation 
was  Increased  from  0.001  to  0.002.  In  addition,  to  account  for  tine -varying, 
random  fluctuations  In  wave  height  that  tend  to  create  relatively  smooth 
macroscale  features  observed  on  the  subaqueous  profile  In  nature,  a  30  percent 
variation  In  recorded  wave  height  was  Imposed  on  the  entire  time  series.  A 
time  step  of  5  min  and  a  length  step  of  2  m  were  used  for  all  simulations. 

164.  Finally,  It  was  Inferred  that  water  level  measured  at  the  gage 
(located  approximately  1  km  offshore)  vas  not  entirely  representative  of  water 
elevations  at  the  shoreline  Influencing  profile  change.  Calculated  beach 
profile  response  best  matched  measured  profile  change  when  water  level  at  the 
gage  was  Increased  by  approximately  0.3  m  during  the  storm  event.  This  Is 
consistent  with  onshore-  and  offshore -directed  wind  stress  associated  with 
storms  affecting  the  southern  California  coast  that  typically  create  fluctua¬ 
tions  In  water  surface  elevation  at  the  shoreline,  producing  setup  as  a  storm 
approaches  and  setdown  as  a  storm  passes. 

Storm  Event  1 

165.  Figure  20  shows  the  time  series  of  wave  and  water  level  character¬ 
istics  for  the  period  29  December  1975  through  9  January  1976.  Measured  wave 
height,  peak  wave  period,  and  water  level  were  available  at  6-hr  Intervals. 
Pressure  data  from  the  array  of  sensors  were  depth-corrected  and  converted  to 
sea-surface  elevation,  yielding  the  significant  wave  height.  Since  the  time 
step  used  In  the  numerical  model  (typically  5  min)  was  shorter  than  the  wave 
and  water  level  measurement  interval,  input  values  for  time  steps  between 
measurements  were  Interpolated. 

166.  A  2-day  period  near  the  beginning  of  the  record  Illustrates  storm 
energy  resulting  from  a  relatively  weak  North  Pacific  low  pressure  system. 
Maximum  wave  heights  were  approximately  1.8  m  with  a  peak  period  of  11.1  sec. 
Elevated  water  levels  due  to  storm  surge,  coincident  with  high  tide,  produced 
a  measured  water  level  at  the  gage  of  1.2  m  at  the  time  of  maximum  wave 
height.  Post-storm  wave  height  decreased  to  less  than  0.9  m,  and  peak  wave 
period  Increased  as  swell  wave  energy  became  dominant. 

167.  The  distribution  of  median  grain  diameter  along  the  profile  Is 
required  to  run  SBEACH.  Grain  size  data  were  not  available  for  this  time 
period  but  median  grain  diameter  recorded  near  the  site  In  October  1983  showed 
little  variation  between  the  1-  and  -1-m  contours  (Corps  of  Engineers  1984). 
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Time  from  0900  12/29/75  (days) 


figure  20.  Tine  series  of  wave  and  water  level  characteristics, 
winter  storm  1975 
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Therefore,  an  average  median  grain  diameter  of  0.17  mm  was  used  for  model 
calibration. 

168  Since  this  storm  was  of  minimal  Intensity,  recorded  beach  profile 
erosion  was  relatively  minor.  Figure  21  illustrates  measured  and  calculated 
profile  shape  for  the  11-day  period.  Calculated  profile  change  from  29 
December  Co  31  December  1973  was  concentrated  on  the  upper  half  of  the 
foreshore  and  was  related  to  locally  elevated  water  levels.  Generally, 
calculated  profile  response  above  NGVD  was  slightly  less  than  measured 
~  adjustments;  the  2-m  contour  retreated  approximately  4.6m  and  profile  volume 
change  above  the  mean  shoreline  (y  -  0)  was  -8  a?/n.  Only  small  adjustments 
In  bar  elevation  were  recorded. 

169.  Measured  profile  response  between  29  December  1973  and  9  January 
1976  Includes  effects  of  the  2-day  storm  event  as  well  as  approximately  8  days 
of  post-storm  recovery  wave  action.  After  31  December  1975,  most  profile 
changes  occurred  below  NGVD  where  two  new  bar  features  developed.  However, 
cumulative  retreat  of  the  2-m  contour  increased  to  10  m,  and  associated 
subaerial  profile  volume  change  increased  to  -IS  a?/m.  The  calculated  profile 
for  9  January  1976  shows  a  small  bar  located  at  about  the  0.6-m  depth  contour, 
In  close  proximity  to  the  Innermost  measured  bar.  Although  the  measured  bar 
Immediately  seaward  of  this  position  Is  not  represented  on  the  calculated 
profile,  the  broad,  low-relief  deposit  at  about  the  -2.4-m  contour  Is  fairly 
well  represented  by  the  simulation.  Calculated  profile  change  above  NGVD 
showed  increased  retreat  of  the  2-m  contour  between  31  December  1975  and 

9  January  1976,  resulting  In  a  cumulative  retreat  of  about  6.4  m.  Simulated 
profile  volume  adjustments  gave  a  cumulative  loss  of  14.8  n?/m.  A  comparison 
of  measured  and  calculated  profile  changes  shows  good  agreement  for  this  storm 
event . 

S&orrc . Event-2 

170.  Verification  of  SBEACH  simulations  was  accomplished  using  process 
and  response  data  from  a  second  small  storm  event  which  occurred  at  Torrey 
Pines  Beach  approximately  1  year  after  Storm  Event  1.  Data  were  analyzed  and 
prepared  for  model  application  In  a  manner  Identical  to  that  for  the  winter 
1975  storm.  Figure  22  shows  the  time  series  of  wave  and  water  level  para¬ 
meters  for  the  period  21  December  1976  to  4  January  1977.  In  this  case,  the 
record  covers  14  days  with  an  increase  in  wave  height  occurring  toward  the 
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Elevation  Relative  to  NGVD  (m)  Elevation  Relative  to  NGVD  (m) 
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Torrey  Pines,  California 


Time  from  1200  12/21/78  (days) 

Figure  22.  Tine  series  of  wave  and  wacer  level  characteristics , 
winter  1976  atom 
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middle  portion  of  the  time  series.  Small  variations  in  wave  period  suggest 
that  locally  generated  winds  likely  were  not  the  main  source  of  energy  for 
this  event.  Consequently,  an  Increase  in  water  level  due  to  wind  setup  at  the 
shoreline  was  not  Imposed  and  only  direct  measurements  of  water  level  at  the 
gage  were  used  for  simulations.  All  other  previously-determined  model 
parameters  were  held  constant. 

171.  A  maximum  wave  height  of  1.6  m  was  recorded  on  27  December  1976 
associated  with  a  15.1-sec  wave  period.  For  the  next  4.5  days  wave  height 
averaged  about  1.1  m.  Fluctuations  in  water  level  were  at  a  minimum  during 
this  time;  however,  water  surface  elevation  remained  above  NGVD  while  wave 
heights  were  greatest.  This  characteristic  was  not  present  in  Storm  Event  1 
and  likely  contributed  to  an  increase  in  eroded  volume  above  NGVD  during  Storm 
Event  2.  Wave  energy  associated  with  North  Pacific  swell  was  dominant  in  the 
14 -day  record. 

172.  This  second  event  exhibited  lower  wave  energy  than  the  first 
storm;  however,  changes  in  profile  volume  and  contour  position  above  NGVD  were 
slightly  greater.  This  can  be  attributed  to  the  duration  of  water  level  above 
the  mean  water  line  during  higher  wave  conditions.  Figure  23  illustrates 
measured  and  calculated  profile  change  for  the  period  21  December  1976  to  4 
January  1977.  Unlike  responses  recorded  during  the  first  storm  event,  a  rela¬ 
tively  thin  veneer  of  sand  was  eroded  from  the  original  profile  betveen  the 
berm  crest  (approximately  2.3-m  contour)  and  the  -2.1-m  contour.  This  volume 
of  sand  was  not  conserved  along  the  profile  since  deposits  of  near-equal 
volume  were  not  recorded  at  a  seaward  position.  Redistribution  of  sand  by 
longshore  processes  is  a  likely  cause  for  the  lack  of  cross -shore  conservation 
of  mass  during  this  time  period  of  dominant  swell  energy.  Measurements  of 
beach  profile  change  above  NGVD  showed  a  maximum  retreat  associated  with  the 
O-depth  contour  (12.8  m  recession),  and  a  steady  decrease  in  contour  retreat 
at  higher  elevations.  The  1-m  contour  receded  8.8  m  and  the  2-m  contour 
receded  4  m.  Measured  changes  in  profile  volume  above  NGVD  were  approximately 
-16.3  mV“- 

173.  Calculated  profile  change  satisfactorily  verified  predictive 
capabilities  of  the  model.  Changes  in  position  of  the  1-  and  2-m  contours 
were  nearly  identical  to  measured  changes  (-10.4  m  and  -4  m) ,  and  the  0-depth 
contour  was  slightly  underpredicted.  Therefore,  profile  volume  change  above 
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NGVD  was  consistent  with  measured  changes  (-17.6  v?/m) .  However,  calculated 
profile  shape  on  the  seawardmost  portion  of  the  active  profile  was  signifi¬ 
cantly  different  than  measured  change.  This  apparent  Inconsistency  is 
explained  by  examining  the  way  in  which  profile  change  is  calculated.  At  the 
shoreward  end  of  the  profile,  the  runup  limit  constitutes  a  boundary  where 
cross-shore  sand  transport  becomes  zero.  The  seaward  boundary  is  determined 
by  the  depth  at  which  no  significant  sediment  motion  occurs.  Once  the  distri¬ 
bution  of  the  sand  transport  rate  is  known,  profile  change  is  calculated  from 
the  mass  conservation  equation.  If  the  assumption  of  conservation  of  mass  in 
a  two-dimensional  sense  is  not  approximated  by  a  three-dimensional  natural 
process,  either  an  excess  or  deficit  of  sand  volume  will  be  predicted. 
Verification  of  calculated  profile  response  for  this  region  vas  otherwise  good. 
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PART  VI:  MODEL  APPLICATION  TO  BEACH  FILL  DESIGN 


174.  Beach  nourishment  Is  becoming  a  preferred  shore  protection  measure 
because  it  is  often  the  least  expensive  solution,  is  flexible  and  can  augment 
existing  shore  protection  structures,  and  provides  a  natural  and  enjoyable 
coastal  environment.  Prediction  of  the  behavior  of  beach  fill  is  required  to 
determine  the  benefits  of  a  nourishment  project  in  relation  to  other  shore 
protection  alternatives  as  well  as  to  estimate  its  longevity  and  function¬ 
ality.  A  realistic  description  of  coastal  sediment  transport  processes  and 

.beach  fill  change  Involves  numerous  hydrodynamic,  sediment,  and  morphologic 
characteristics  which  vary  in  space  and  time.  Numerical  modeling  of  beach 
profile  and  fill  change  provides  a  general  and  quantitative  means  of  repre¬ 
senting  the  combined  and  simultaneous  effects  of  these  variables.  Application 
to  beach  fill  design  was  one  of  the  objectives  in  developing  SBEACH. 

175.  Figure  24  is  a  definition  sketch  of  a  nearshore  region  containing 
a  beach  fill.  The  fill  is  acted  on  by  ordinary  cross-shore  wave  processes, 
such  as  those  associated  with  episodic  and  potentially  catastrophic  storms, 
which  mold  it  from  an  Initial  shape  to  a  dynamic  equilibrium  shape.  The  fill 
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Figure  24.  Definition  sketch  categorizing  a  beach  fill 
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is  also  exposed  to  longshore  currents  induced  by  waves  and  wind,  which  move 
sand  alongshore  and  deplete  the  fill  when  updrift  sand  source  is  available  or 
maintenance  nourishment  is  not  provided.  Longshore  transport  does  not  enter 
in  SBEACH;  therefore,  the  model  is  most  appropriately  applied  to  the  internal 
region  of  a  long  fill  and  far  from  coastal  structures  and  boundaries  where 
differentials  in  longshore  transport  may  be  large. 

176.  This  chapter  describes  hypothetical  applications  of  SBEACH  to 
predict  the  initial  adjustment  of  beach  fill  to  wave  action  by  storms. 

-  Initially,  simulations  are  made  for  two  types  of  storms  occurring  on  the  US 
East  Coast.  One  synthetic  tropical  storm  (hurricane)  and  one  synthetic 
extratropical  storm  (northeaster)  reoresenting  mid-Atlantic  events  with 
approximate  2  to  5-year  return  periods  are  used  to  predict  erosion  associated 
with  two  hypothetical  beach  fill  configurations  and  subsequent  post-storm 
■ecovery.  Eroded  sand  volume  and  contour  movement  are  evaluated  as  a  function 
of  storm  type,  fill  cross  section,  grain  size,  and  time. 

177.  In  the  second  half  of  the  chapter,  simulations  are  made  of  the 
response  of  two  different  fill  cross  sections  of  two  different  grain  sizes  to 
a  tropical  storm  with  a  return  period  of  about  10  years.  Comparisons  are 
made  to  a  beach  without  a  fill  to  evaluate  the  effectiveness  of  artificial 
nourishment  as  a  means  of  shore  protection.  A  seawall  is  included  in  the 
simulations  to  illustrate  model  capabilities  to  predict  profile  evolution  in 
situations  involving  a  vertical  wall.  Post-storm  recovery  is  also  simulated 
for  this  case  making  it  feasible  to  assess  the  possibility  for  the  eroded  fill 
sand  to  contribute  material  in  the  beach  recovery  process.  Some  general 
guidelines  are  provided  concerning  the  design  of  beach  nourishment  schemes 
based  on  the  numerical  simulations. 

178.  Finally,  a  concluding  discussion  is  given  regarding  capabilities 
and  limitations  of  the  numerical  model  as  a  tool  for  predicting  beach  fill 
evolution.  Much  of  the  r.aterial  presented  in  this  chapter  may  be  found  in 
Kraus  and  Larson  (1989)  ano  Larson  and  Kraus  (1989b,  in  press). 

Fill  Response  to  Moderate  Storms 

179.  A  verified  model  enables  examination  of  the  shore  protection 
functioning  of  alternative  fill  cross  sections  or  fill  templates  at  a  beach 
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nourishment  project  impacted  by  storms.  Such  calculations  must  be  made  to 
perform  cost-benefit  analysis  and  to  estimate  the  minimum  fill  volume  needed 
to  provide  a  specified  level  of  storm  protection  against  erosion  and  inunda¬ 
tion.  In  an  overall  evaluation,  longshore  variability  of  the  protective  beach 
complex  as  well  as  storms  must  be  considered  (van  de  Graaff  1983,  Savage  and 
Blrkemeier  1987) . 

180.  Two  approaches  can  be  taken  to  estimate  storm  impact  on  a  beach 
fill  complex;  one  may  be  called  the  design- storm  approach  and  the  other  the 
storm-ensemble  approach.  The  design  storm  is  either  a  hypothetical  or  a 
.historical  event  which  produces  a  specified  storm  surge  hydrograph  and  wave 
condition  at  the  project  site.  Surge  is  a  water  level  rise  caused  by  wind 
stress  and  atmospheric  pressure  variation;  waves  also  produce  a  rise  in  mean 
water  level  through  setup.  The  time  average  (on  the  order  of  an  hour)  of 
surge,  wave  setup,  and  tide  is  called  the  stage.  In  a  stage -frequency 
analysis,  the  design  storm  may  have  a  certain  frequency  of  occurrence  or 
return  period  (for  example,  a  100-year  storm). 

181.  The  design  storm  approach  is  problematic  fo>  „se  in  beach  erosion 
modeling  because  profile  change  is  sensitive  to  storm  duration,  surge  shape, 
and  wave  height  and  period,  in  addition  to  peak  stage.  The  maximum  water 
level  associated  with  the  surge  of  a  design  storm  may  produce  less  erosion 
than  a  storm  of  lower  surge  but  longer  duration,  or  a  storm  of  lower  surge  but 
higher  waves.  However,  in  the  absence  of  an  extensive  storm  data  base,  use  of 
several  large  historical  storms  (large  in  the  sense  of  having  caused  signifi¬ 
cant  erosion)  to  drive  the  model  is  considered  acceptable.  In  this  approach, 
it  is  recommended  that  the  phasing  of  the  peak  tide  be  varied  with  the  peak 
surge  since  model  results  are  sensitive  to  water  level.  For  example,  arrival 
of  the  peak  surge  during  low  tide  on  a  mesotidal  or  macrotldal  coast  may 
produce  little  or  no  erosion,  whereas  arrival  of  the  peak  surge  at  high  tide 
could  bring  disaster. 

182.  The  most  satisfactory  solution  to  the  problem  of  the  many-to-one 
relation  between  beach  erosion  and  stage  frequency  is  to  use  the  storm- 
ensemble  approach,  i.e.,  calculate  erosion  for  a  large  number  of  storms  and 
key  the  erosion  to  the  frequency  of  storm  occurrence  (Kriebel  1982;  Scheffner 
1988,  1989).  This  method  yields  an  erosion  or  recession  frequency  of  occur¬ 
rence  curve.  In  this  method,  tropical  and  extratropical  storms  must  be 
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treated  independently  because  they  have  different  physical  characteristics. 
Tropical  storms  are  infrequent  events  of  short  duration  and  high  intensity, 
whereas  extratropical  storms  are  more  frequent  and  usually  of  longer  duration 
and  lower  intensity.  The  storm-ensemble  approach  is  recommended  for  project 
design,  although  it  requires  a  storm  data  base  and  is  much  more  computation¬ 
ally  intensive  than  the  design-storm  approach. 

183.  In  this  part  of  the  report,  the  response  of  nourished  beaches  to  a 
representative  tropical  storm  and  an  extratropical  storm  is  calculated  to 
examine  predictions  of  the  model  to  storms  of  differing  waves  and  water  levels 
(surge  and  duration).  The  two  storms  were  synthesized  to  produce  erosion 
resulting  from  a  2  to  5  year  storm  surge  for  the  mid-Atlantic  Ocean  coast  of 
the  United  States.  The  amount  of  eroded  volume  for  such  events  Is  on  the 
order  of  20  to  30  m^/m  (Savage  and  Blrkemeier  1987). 

Design  storms 

184.  Surge  hydrographs  developed  for  the  hypothetical  storms  are  shown 
in  Figure  25.  The  hurricane  surge  has  a  duration  of  approximately  12  hr,  with 
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a  peak  surge  of  2  a,  and  a  duration  above  half  the  peak  surge  (1  m)  of  6  hr. 
The  shape  of  the  hydrograph  follows  an  inverse  hyperbolic  cosine  squared.  The 
surge  of  the  northeaster  has  a  duration  of  36  hr,  with  a  peak  surge  of  1  a, 
and  a  duration  above  half  the  peak  surge  (0.5  m)  of  18  hr.  A  cosine-squared 
function  was  used  to  generate  the  shape  of  the  northeaster  surge.  The  peak 
surge  of  the  hurricane  is  higher  because  wind  speeds  in  hurricanes  are  on 
average  greater  than  those  in  northeasters. 

185.  The  tiae  histories  of  wave  height  and  period  assigned  to  the 
storms  are  shown  in  Figure  26.  Both  have  peak  deep  water  wave  heights  of  5  m, 
which  occur  during  the  time  when  the  respective  surges  are  greater  than  half 
the  maximum.  The  duration  of  high  waves  for  the  northeaster  is  thus  three 
times  that  of  the  hurricane.  Since  the  radius  of  a  northeaster  is  typically 
several  times  greater  than  that  of  a  hurricane,  the  fetch  is  longer,  resulting 
in  longer  wave  periods.  Wave  height  and  period  of  1  m  and  7  sec  were  applied 
for  6.5  days  before  the  start  of  the  storms  to  mold  the  profiles  into  a 
realistic  shape  under  ordinary  waves.  As  the  storm  passed,  wave  height  and 
period  were  changed  to  0.5  m  and  10  sec  to  simulate  long-period  recovery  swell 
wave  conditions.  A  sinusoidal  tide  was  applied  with  a  12-hr  period  and  0.5-m 
amplitude;  a  peak  in  the  tide  occurred  during  the  peak  surge  of  each  storm. 
Mach  fill  shapes 

186.  Following  the  procedure  of  Kraus  and  Larson  (1989),  two  different 
beach  fill  cross  sections  or  templates  were  designed  for  exposure  to  storms. 
The  artificial  berm  had  most  of  the  fill  placed  on  the  beach  and  above  mean 
sea  level.  The  as -constructed  berm  extended  horizontally  for  a  distance  of 
16  m  at  an  elevation  of  3  m  and  then  tapered  with  a  1:20  slope  to  join  the 
original  beach  profile  at  a  depth  of  1.4  m  (Figure  27).  This  is  a  common 
beach  fill  design  in  the  United  States  (e.g.,  see  Jarrett  and  Hemsley  1988). 
The  ocher  fill  template  is  termed  profile  nourishment,  a  concept  advocated  by 
Bruun  (1988).  In  the  example,  the  material  was  placed  over  the  profile  from 
an  elevation  of  1  m  to  -2  m  (Figure  28)  in  an  approximation  of  the  existing 
profile  shape.  The  fill  volume  was  the  same  for  each  template  (140  m^/a)  . 

187.  Initially,  all  simulations  were  made  assuming  a  0.2-mm  sand  beach, 
with  the  fills  and  the  natural  beach  having  the  same  grain  size.  However, 
simulations  were  also  made  for  fill  grain  sizes  in  0.1-mm  increments  from  0.2 
to  1.0  mm.  In  these  cases,  the  grain  size  was  specified  as  that  of  the  fill 
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Figure  26.  Time  history  of  wave  height  and  period  for  the  synthetic  storms 
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Figure  27.  Response  of  artificial  bens  to  hurricane  impact 
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over  the  portion  of  the  profile  originally  occupied  by  the  fill,  and  0.2  ma 
elsewhere.  Profile  elevations  are  given  with  respect  Co  mean  sea  level,  and  a 
water  temperature  of  20  deg  was  specified  in  the  model  for  computation  of  the 
sand  fall  speed. 

Profile  chance 

188.  Figures  27  through  30  illustrate  predicted  impacts  of  the  two 
storms  on  the  two  nourishment  schemes.  The  curve  labelled  "profile  without 
fill’  gives  a  hypothetical  dune,  beach,  and  subaqueous  equilibrium  profile  for 

"reference.  The  construction  profile  configuration  is  denoted  as  "profile  with 
fill,"  whereas  the  "pre- storm  profile"  shows  the  profile  after  6.5  days  of 
ordinary  waves.  The  "post- storm  profile*  represents  the  shape  prior  to  the 
recovery  wave  period  and  the  "post-storm  recovery"  profile  displays  the 
profile  after  experiencing  approximately  two  weeks  of  recovery  waves.  The 
following  few  paragraphs  discuss  pre-  and  post-storm  profile  development. 

189.  Pre -storm.  Pre-storm  profiles  ("design  profile"  of  Figure  24)  of 
the  artificial  berm  and  nourished  profile  differ  significantly  in  the  inner 
surf  zone.  A  pronounced  scarp  is  produced  with  the  artificial  berm,  whereas 
the  nourished  profile  experiences  more  subtle  changes  since  it  was  placed  in  a 
near-equilibrium  configuration.  For  both  cross  sections,  a  small  breakpoint 
bar  formed  at  about  the  210-  to  220-m  mark  (measured  from  an  arbitrary 
baseline).  Sand  was  removed  from  the  foreshore  and  distributed  along  the 
profile  just  beyond  the  2-m  depth  contour.  Thus,  regardless  of  the  initial 
fill  configuration,  the  model  predicts  that  fill  material  will  be  transported 
a  considerable  distance  offshore  to  the  point  of  incipient  wave  breaking  in 
the  process  of  molding  the  surf  zone  profile  to  an  equilibrium  shape. 

190.  Post-storm.  Subaqueous  portions  of  the  post-storm  profiles  are 
quite  similar  in  shape,  being  reworked  by  storm  waves  to  the  same  equilibrium 
shape  independent  of  initial  profile  configuration  and  type  of  storm.  The 
shoreline  position  (0-depth  contour)  actually  moved  seaward  of  its  pre-storm 
location,  with  the  accreted  material  supplied  from  the  ordinarily  subaerial 
portion  of  the  profile  that  was  Inundated  due  to  storm  surge.  A  small  bar 
formed  at  approximately  the  5 -a  depth  contour  under  the  high  storm  waves,  but 
it  is  not  shown  here  to  better  display  changes  near  the  beach. 

191.  An  Important  outcome  of  the  predictions  is  that,  under  the  action 
of  simulated  hurricane  and  northeaster  events,  resultant  profile  change  was 
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Figure  29.  Response  of  artificial  berm  to  northeaster  impact 
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Figure  30.  Response  of  nourished  profile  to  northeaster  Impact 
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very  similar.  The  use  of  one  typical  storm  descriptor  (e.g.  maximum  stage)  to 
estimate  shoreline  recession  or  eroded  volume  would  have  predicted  greater 
changes  associated  with  the  hurricane,  thus  producing  misleading  results. 

192.  Post-storm  recovery.  In  all  cases,  substantial  accretion  occurred 
near  MSL  during  the  post-storm  recovery  phase.  The  upper  foreshore  of  the 
nourished  profile  experienced  more  accretion  than  the  artificial  berm  cases. 
These  results  are  consistent  with  the  concept  that  a  beach  profile  with  a 
natural  shape  can  best  respond  to  changes  In  Incident  waves. 

~Eroded  volume  and  contour  change 

193.  Figure  31  shows  the  evolution  of  eroded  volume  above  the  0-  and 
1-m  contours  (and  0.5-m  contour  for  the  profile  nourishment  case).  Eroded 
volume  above  the  0-depth  contour  Increases  rapidly  at  the  beginning  of  pre¬ 
storm  ("ordinary")  wave  action,  describing  the  behavior  of  fill  material 
during  initial  profile  adjustment.  In  contrast,  eroded  volume  above  the  1-m 
contour  shows  a  much  slower  increase.  The  profile  nourishment  case  Illu¬ 
strates  greater  initial  erosion  during  the  early  stage  of  wave  action,  but 
also  greater  recovery  during  the  post-storm  period.  It  should  be  emphasized 
that  longshore  processes  are  omitted  from  the  present  discussion;  however, 
beach  nourishment  places  a  greater  amount  of  material  Into  the  active  littoral 
zone,  thereby  Increasing  its  potential  for  transport  alongshore  and  out  of  the 
project  reach. 

194.  The  volume  of  eroded  material  above  the  0-depth  contour  does  not 
show  a  significant  increase  during  the  storms,  changing  from  about  22  to 

27  m^/m  for  erosion  of  the  artificial  berm  during  the  northeaster.  Eroded 
volume  above  the  1-m  contour  abruptly  increases  at  the  start  of  the  storms, 
from  about  3  to  17  m^/m  for  the  artificial  berm  Impacted  by  the  northeaster. 
The  volume  of  erosion  above  the  0-depth  contour  Is  relatively  unchanged 
because  these  moderate  storms  primarily  remove  material  from  the  upper  portion 
of  the  profile  and  redistribute  it  over  the  shallow  nearshore  region. 

195.  The  hurricane  and  northeaster  produce  nearly  the  same  amount  of 
erosion,  25  to  30  m^/  above  the  0-depth  contour  and  13  to  16  m^/m  above  the 
1-m  contour.  Therei.„..e,  the  longer  surge  duration  of  the  northeaster  was 
approximately  equivalent  In  erosion  capacity  to  the  higher  surge  of  the 
shorter  duration  hurricane.  Eroded  volumes  above  the  0-depth  contour  are 
comparable  to  those  associated  with  2  to  5  year  return  period  storms  impacting 
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Figure  31.  Eroded  volume  above  specific  ctntours 
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the  mid-Atlantic  coast  (Savage  and  Birkemeier  1987).  Evolution  of  the  eroded 
volumes  above  mean  sea  level  shows  an  approximate  exponential  approach  to  an 
equilibrium  value.  These  trends  are  In  general  agreement  with  those  obtained 
by  Krlebel  and  Dean  (1985),  who  numerically  examined  eroded  volume  and  berm 
recession  as  a  function  of  wave  height,  surge  level,  and  other  parameters. 

196.  It  Is  Interesting  that  there  is  a  tendency  for  the  hurricane- 
impacted  beach  to  recover  more  rapidly  than  the  northeaster- impacted  beach 
(Figure  31).  This  is  In  agreement  with  qualitative  observations  of  post-storm 
recovery  along  the  Florida  coast.  An  explanation  lies  in  the  nature  of  the 
storm  hydrograph;  material  Is  not  removed  as  far  offshore  during  hurricanes  as 
in  the  case  of  the  longer  duration  and  lower  surge  extratropical  storms.  This 
material  is  deposited  closer  to  shore  and  can  be  moved  back  on  the  beach  more 
rapidly  under  recovery  wave  action. 

197.  Figure  32  displays  the  evolution  of  the  position  of  the  0-  and  1-m 
contours.  A  decrease  in  contour  position  indicates  recession  of  the  beach  at 
that  contour.  The  0-depth  contour  for  both  fill  types  and  both  storm  types 
shows  recession  during  pre-storm  and  storm  periods,  but  begins  to  advance 
before  the  end  of  the  storms  and  prior  to  arrival  of  recovery  waves,  as  the 
surge  subsides  and  the  wave  height  decreases.  The  1-m  contour  shows  no 
recovery  for  the  artificial  berm  because  the  foreshore  is  relatively  steep, 
whereas  the  0.5-m  contour  for  the  nourished  profile  does  show  some  recovery 
since  the  gentler  slope  allows  post-storm  wave  runup  to  build  a  berm.  In 
SBEACH,  berm  formation  and  growth  are  largely  controlled  by  the  elevation 
reached  by  wave  runup  (Larson  1988,  Report  1),  and  accretionary  processes  are 
diminished  on  steeper  slopes. 

198.  Figure  33  Illustrates  eroded  volume  at  the  end  of  the  storm  (prior 
to  recovery  wave  action)  as  a  function  of  beach  fill  grain  sire.  As  previously 
mentioned,  median  grain  diameter  of  the  native  beach  was  set  at  0.2  mm,  and  the 
area  in  which  the  fill  was  placed  was  assigned  the  grain  size  of  the  fill. 

This  procedure  does  not  allow  tracking  of  movement  of  the  different  grain 
sizes.  However,  since  surf  zone  sediment  is  usually  sorted  with  coarser 
material  located  higher  on  the  active  profile,  this  simple  procedure  is 
considered  to  provide  a  reasonable  first  approximation  of  the  response  of  a 
natural  beach  of  varying  grain  size. 
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199.  For  the  particular  synthetic  storms  used,  Figure  33  shows  a  rela¬ 
tively  steep  decrease  in  eroded  volume  as  grain  size  increases  from  0.2  to 
0.4  mm,  with  a  gentle  decrease  thereafter.  This  behavior  results  from  the 
property  of  the  empirically  determined  functional  dependence  of  wave  energy 
dissipation  needed  to  generate  an  equilibrium  profile  of  a  given  grain  size. 
The  property  Is  manifested  by  computed  dissipation  rates  which  rise  steeply  in 
the  range  of  0.1  to  0.4  mm,  and  then  increase  at  a  lower  rate  with  Increasing 
grain  size  (Moore  1982).  Since  the  rate  of  decrease  In  erosion  Is  small 
beyond  0.4  mm,  and  the  cost  of  beach  fill  typically  Increases  substantially 
for  larger  size  material,  calculations  such  as  those  illustrated  In  Figure  33 
allow  an  evaluation  to  be  made  of  initial  nourishment  volume  and  subsequent 
fill  maintenance  costs.  In  project  planning,  curves  similar  to  those  in 
Figure  33  should  be  developed  with  consideration  of  the  design  profile  and  the 
storm  characteristics  of  the  area,  using  the  storm-ensemble  approach. 

Beach  Fill  Response  With  a  Seawall 

200.  Another  aspect  of  the  simulation  was  to  display  the  capabilities 
of  the  model  to  predict  profile  evolution  of  a  fill  fronting  a  seawall.  A 
schematic  storm  event  was  developed  similar  to  the  previous  two  synthetic 
storms  and  consisted  of  a  period  of  low  waves  when  profile  adjustment  took 
place,  a  3-day  storm  event,  and  a  period  of  recovery  waves.  The  total 
simulation  time  was  30  days  and  comparisons  were  made  for  a  beach  with  and 
without  fill.  The  two  fill  types  studied  were  again  an  artificial  berm  and 
profile  nourishment,  with  cases  examined  using  median  grain  diameters  of  0.25 
and  0.40  mm. 

201.  The  first  21  days  rf  the  simulation  consisted  of  constant  waves 
and  water  level,  with  a  wave  height  of  0.5  m  and  a  wave  period  of  8  sec. 

These  values  were  assumed  to  represent  typical  wave  conditions  on  an  open- 
ocean  Atlantic  coast  ind  were  used  to  produce  a  realistic  equilibrium  profile. 
Under  these  conditions,  beach  change  can  be  mildly  eroslonal  or  accretionary, 
depending  on  the  grain  size  and  Initial  profile  shape.  At  day  21  ("Pre-storm" 
in  Figure  34)  a  storm  was  imposed  that  lasted  for  3  days.  The  average  surge 
water  level  during  the  storm  followed  a  Gaussian  form  and  reached  a  maximum 
elevation  of  3  m  above  MSL.  The  average  wave  height  was  given  a  sinusoidal 
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Figure  34.  Time  history  of  waves  and  water  level 


variation  during  the  storm,  having  a  peak  of  2.5  m  at  the  peak  of  the  surge 
( "Mid- storm" ) .  Wave  period  during  the  3-day  event  varied  between  8  and  12 
sec,  with  12-sec  waves  occurring  at  the  time  of  the  maximum  wave  height.  At 
the  end  of  the  storm  ("Post-storm"),  the  surge  was  considered  non-existent  and 
waves  arrived  as  swell  with  a  height  of  0.5  m  and  a  period  of  16  sec. 

202.  In  the  following  sections,  adjustment  of  the  profile  to  the  30-day 
design  sea  conditions  is  examined  with  SBEACH  for  three  situations:  an 
original  (existing)  condition,  a  fill  of  artificial  berm-type,  and  a  fill 
carried  out  using  the  profile  nourishment  technique.  The  simulations  were 
done  using  parameters  from  the  calibrated  model  with  a  time  step  of  20  min  and 
a  space  step  of  5  m. 

Existing  conditions 

203.  The  profile  configuration  representing  the  beach  prior  to  nourish¬ 
ment  activity  is  shown  in  Figure  35.  A  2-m  high  foreshore  is  backed  by  an 
infinitely  high  seawall,  although  its  crown  is  truncated  at  a  2-m  elevation  in 
this  and  subsequent  figures.  The  original  profile  decreases  linearly  to  the 
O.S-m  depth  contour  with  two  different  slopes;  the  break  in  slope  located  at 
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Figure  35.  Simulation  of  profile  evolution  for  existing  beach 
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an  elevation  of  1  m.  Seaward  of  the  segmented  linear  profile,  an  equilibrium 
shape  is  used  Into  deeper  water. 

204.  Simulations  showed  that  the  0.25-mra  sand  beach  eroded  on  the 
foreshore  during  the  first  21  days  of  low  waves,  and  a  small  bar  formed  with  a 
crest  located  about  120  m  from  the  base  line.  In  contrast,  the  0.40-mm  beach 
accreted  slightly  on  the  foreshore,  indicating  chat  the  transport  was  directed 
onshore.  At  mid* storm,  two  bars  formed  offshore,  the  berm  was  completely 
submerged,  and  the  beach  eroded  considerably  at  the  seawall  for  both  sand 
sizes.  The  inner  bar  was  created  at  the  beginning  of  the  surge,  whereas 
formation  of  the  outer  bar  corresponded  to  the  peak  surge  and  maximum  wave 
height.  Post-storm  profile  shape  shoved  that  the  berm  had  eroded  further  at 
the  seawall,  slightly  more  so  for  the  0.25-mm  sand  beach  than  for  the  0.40-mm 
sand  beach,  and  a  single  large  bar  appeared  on  both  profiles.  The  bar  on  the 
finer  sand  beach  is  broader  and  lower  in  elevation  with  respect  to  the  initial 
profile  than  the  coarser  sand  bar.  It  is  interesting  to  note  that  the 
position  of  the  post-storm  shoreline  is  at  about  the  same  location  (approx¬ 
imately  75  m  from  baseline)  for  both  beaches;  however,  sand  moved  further 
offshore  on  the  finer  sand  beach,  and  its  subaerial  section  is  more  deflated. 

205.  Lower-steepness  waves  arriving  during  the  7-day  recovery  period 
passed  over  the  storm  bars  on  both  beaches,  causing  little  change  in  bar 
shape.  During  the  recovery  phase,  swell  broke  in  shallower  water,  creating  a 
small  secondary  bar  Inshore  on  both  profiles.  The  inner  bar  was  formed  by 
onshore  transport,  producing  a  deep  trough  between  the  storm  bar  and  inshore 
bar.  Both  beaches  accreted  to  the  limit  of  runup. 

berm  fill 

206.  Simulation  results  for  the  artificial  3-m  high  berm  are  given  in 
Figure  36.  A  beach  fill  of  volume  85  mV»  with  the  same  grain  size  as  the 
existing  beach  was  placed  mainly  on  the  subaerial  portion  of  the  profile.  The 
fill  extended  horizontally  from  the  seavall  for  20  m,  tapered  down  for  several 
meters,  then  extended  offshore  with  slope  1:10  to  meet  the  equilibrium  profile 
at  a  depth  of  approximately  1  m. 

207.  Twenty-one  days  of  pre-storm  waves  eroded  the  toe  of  the  0.25-mm 
artificial  sand  berm  and  created  a  small  bar,  whereas  the  0.40-mm  beach 
configuration  was  effectively  stable  under  the  mild  waves.  Berm  elevation  for 
both  sand  sizes  was  lowered  during  the  storm,  but  not  to  the  extent  of  the 
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original  beach  case  (Figure  35).  The  artificial  berm  provided  substantial 
protection,  but  should  have  been  higher  to  prevent  waves  from  reaching  the 
seawall  for  this  extreme  storm  surge.  The  pre-storm  bar  was  located  closer  to 
shore  than  in  the  existing  case  simulation.  Post-storm  bars  for  both  sand 

-  sizes  were  located  in  essentially  the  same  position  as  for  the  corresponding 
existing  beach  cases. 

208.  The  recovery  cross  sections  for  both  sand  sizes  were  considerably 
different  than  the  original  beach  cases  because  of  the  greater  amount  of  sand 

-in  the  system.  Although  significant  recovery  took  place  for  the  0.2S-ma  sand 
beach,  sand  was  effectively  trapped  in  the  storm  bar  and  could  not  contribute 
to  post-storm  beach  recovery.  Eroded  material  for  the  coarser  sand  beach  was 
deposited  closer  to  shore  and,  during  the  recovery  phase,  the  entire  storm 
profile  translated  shoreward.  Although  a  limited  amount  of  foreshore  recovery 
took  place  with  the  0.40-mit  sand,  a  significant  volume  of  material  remained 
near  to  shore. 

Profile  nourishment 

209.  The  profile  nourishment  case  is  depicted  in  Figure  37  where  the 
major  amount  of  material  is  placed  along  the  nearshore  portion  of  the  profile 
in  accordance  with  the  concept  that  a  beach  can  best  protect  Itself  if  in  a 
natural  shape.  The  same  volume  of  fill  as  in  the  artificial  berm  case 

(85  m^/m)  was  placed  over  the  profile  in  an  equilibrium  shape,  and  tapered  at 
the  seaward  end.  This  fill  template  resulted  in  positioning  of  the  shoreline 
approximately  10  m  seaward  of  the  artificial  berm  case. 

210.  Simulation  results  show  that  subaerial  erosion  was  reduced 
significantly  as  compared  to  existing  conditions,  even  though  the  entire 
profile  was  submerged  for  much  of  the  storm.  Two  bars  developed  again,  but 
the  inner  bar  ("Mid-storm")  for  the  0.25-mm  sand  generated  during  the  rising 
surge  was  located  shoreward  of  the  pre-storm  bar.  Bar  development  (with 
respect  to  the  initial  profile)  for  this  grain  size  was  more  pronounced  than 
that  for  the  artificial  berm  case,  whereas  profile  evolution  with  0.40-mm  sand 
showed  less  bar  development.  The  profile  nourishment  template  showed  greater 
retention  of  sand  above  MSL  than  the  artificial  berm  case,  and  less  material 
was  redistributed  along  the  profile.  However,  because  of  its  low  elevation, 
the  profile  nourishment  template  exposed  the  wall  to  more  wave  action  during 
the  surge  than  did  the  artificial  berm  of  the  same  volume. 
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Figure  37.  Profile  nourishment  simulations 
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nfscusslon  of  Beach  Fill  Simulations 


211.  The  numerical  model  allowed  realistic  time-dependent  calculations 
to  be  made  of  beach  fill  change  in  terms  of  readily  available  engineering 
data.  Model  performance  was  realistic  regarding  storm-induced  beach  erosion, 
showing  correct  dependencies  on  surge  and  wave  input,  as  well  as  initial 
profile  shape.  Storm- Induced  beach  and  dune  erosion  cannot  be  uniquely 
specified  through  a  single  storm-related  parameter  such  as  maximum  stage. 

This  result  demonstrates  the  limited  usefulness  of  the  design  storm  approach. 
Qualitative  features  of  the  post-storm  recovery  phase  were  also  reproduced, 
but  further  work  is  required  to  accurately  simulate  recovery. 

212.  The  model  can  be  used  to  Judge  the  relative  behavior  of  various 
beach  fill  cross  sections  exposed  to  ordinary  and  extreme  waves  for  time 
Intervals  on  the  order  of  days  to  weeks.  Fill  placed  on  the  upper  beach  was 
calculated  to  move  offshore  to  the  region  of  incipient  wave  breaking,  in 
general  agreement  with  the  initial  movement  of  beach  fill  material  (Hansen  in 
preparation).  Also,  a  limited  number  of  calculations  indicated  that  it  may 
not  be  cost-effective  to  use  beach  fill  for  storm  protection  with  a  median 
grain  diameter  much  greater  than  0.4  mm,  owing  to  the  typically  greatly 
Increased  cost  and  the  declining  benefit  in  decreased  volume  of  eroded 
material.  Design  alternatives  for  specific  situations  (fill  grain  size,  fill 
cross  section,  and  storm  climate)  could  be  evaluated  with  the  model. 

213.  The  model  is  based  on  relationships  for  cross-shore  sand  transport 
produced  by  breaking  waves,  implying  that  longshore  sand  transport  is  con¬ 
sidered  negligible.  Thus,  caution  must  be  exercised  In  applying  the  model  and 
Interpreting  results  if  gradients  in  longshore  transport  exist.  Longshore 
transport  is  expected  to  be  a  major  contributing  process  to  beach  change  over 
long  time  periods,  for  example,  seasons  to  a  year.  For  such  time  frames,  the 
present  model  should  be  combined  with  other  predictive  methodologies,  such  as 
a  shoreline  change  model  (e.g.  Kraus  1983;  Kraus,  Hanson,  and  Harlkal  1985; 
Hanson  1989;  Kraus  1989),  which  compute  the  time  rate  of  change  of  shoreline 
position  based  on  estimates  of  wave-induced  longshore  sand  transport. 
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PART  VII;  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 


214.  The  numerical  model  of  beach  profile  change  SBEACH  was  developed 
using  large  wave  tank  data  and  verified  with  field  data,  as  reported  in 
previous  work  (Larson  1988,  Larson  and  Kraus  1989a).  In  the  present  report, 
additional  model  verification  was  made  using  field  data  on  berm  and  dune 
erosion  from  the  US  East  and  Vest  coasts.  Overall  model  performance  for 
simulating  the  erosional  phase  of  a  storm  event  was  satisfactory  whereas  the 
recovery  phase  following  a  storm  event  was  less  well  described  and  essentially 
only  in  qualitative  agreement  with  the  measured  profile  evolution. 

215.  A  high-quality  data  set  describing  storm- induced  erosion  off  the 
coast  of  New  Jersey  was  used  to  Investigate  the  capability  of  the  model  to 
simulate  berm  and  dune  erosion.  SBEACH  predicted  eroded  volumes  and  contour 
retreat  in  agreement  with  what  could  be  confidently  inferred  to  occur  during 
the  erosional  part  of  the  storm  event.  When  the  storm  waned  and  waves  became 
accretionary,  a  large  sand  ridge  formed  on  the  foreshore  of  a  size  that  was 
not  predicted  by  the  model.  However,  qualitative  features  of  the  recovery 
were  reproduced  by  the  model,  including  development  of  profiles  with  different 
shapes  according  to  initial  conditions. 

216.  The  capability  of  the  model  to  describe  profile  change  on  the  US 
West  coast  was  evaluated  with  data  available  from  minor  storms  at  Torrey  Pines 
Beach,  California.  Contour  retreat  on  the  foreshore  was  well  simulated  along 
the  general  profile  shape.  The  amount  of  erosion  predicted  by  the  model  was 
also  in  good  agreement  with  field  measurements.  Similar  quantities  of  calcu¬ 
lated  and  measured  erosion  were  obtained  by  increasing  the  water  level  by 

0.3  m  over  that  measured  at  the  gage  located  offshore  in  9.1  m  of  water.  This 
increase  can  be  ascribed  to  a  number  of  factors  such  as  locally  generated  wind 
setup  landward  of  the  gage,  wave  setup,  and  differences  in  atmospheric 
pressure.  Further  theoretical,  laboratory,  and  field  research  is  needed  to 
better  quantify  hydrodynamic  and  sand  transport  processes  in  the  swash  zone, 
as  model  predictions  are  sensitive  to  changes  in  water  level  nearshore. 

217.  SBEACH  was  applied  to  simulate  the  short-term  fate  of  beach  fill  by 
examining  initial  adjustment  of  the  fill  and  its  response  to  storms.  In  these 
hypothetical  cases,  beach  fill  cross  section  and  grain  size  were  varied,  and 
the  response  of  the  fill  to  synthetic  storm  events  was  investigated.  Synthetic 


103 


storms  corresponding  to  tropical  (hurricane)  and  extratroplcal  (northeasters) 
events  were  chosen  as  being  typical  for  the  US  East  coast.  Two  fill  cross 
sections  were  evaluated,  the  artificial  berm  and  profile  nourishment. 

218.  Simulations  of  fill  response  to  storm  events  showed  the  Importance 
of  both  peak  surge  level,  surge  duration,  and  wave  conditions.  For  example,  a 
high  surge  with  a  low  duration  produced  as  much  erosion  as  a  low-surge  storm 
with  a  longer  duration.  The  profile  nourishment  fill  design  showed  less 
change  for  the  design  storm  than  the  artificial  berm  of  the  same  volume, 
suggesting  that  the  nourished  profile  was  more  In  accordance  with  the  natural 
equilibrium  profile  and  minimized  profile  change  as  the  wave  and  water  level 
conditions  changed.  This  assumption  was  further  substantiated  by  the  stronger 
recovery  shown  by  the  profile-nourished  beach.  Model  simulation  also  indi¬ 
cated  that  little  gain  In  fill  stability  (erosion  protection)  was  achieved 
during  a  storm  if  the  fill  grain  size  exceeded  0.40  mm.  Overall,  the  numeri¬ 
cal  model  provided  reasonable  results  for  simulating  the  response  of  beach 
fills  and  could  be  used  in  a  cost-benefit  analysis  for  determining  the  optimal 
fill  design  and  protection  level.  However,  It  is  necessary  to  verify  the 
model  prior  to  use  in  planning  and  design. 

219.  Conditions  under  which  the  numerical  model  is  expected  to  give 
reliable  results  may  be  summarized  as  follows: 

ft.  Breaking  waves  are  the  main  cause  of  sand  transport  and  thus 
profile  change. 

fe.  Longshore  effects  are  negligible,  that  is,  the  incident  waves 
approach  normal  to  the  shoreline  or  differentials  in  longshore 
transport  are  minor  (influence  of  structures  blocking  longshore 
transport  is  small  and  the  shoreline  is  straight) . 

£.  Median  grain  diameter  on  the  profile  is  in  the  interval  0.1  to 
2.0  mm  and  reasonably  uniform  across-shore. 

220.  Experience  in  simulating  field  situations  has  identified  several 
areas  requiring  basic  research  and  model  refinement.  The  most  important  area 
for  further  research  is  the  description  of  onshore  sand  transport  and  berm 
development.  The  post-storm  recovery  phase  must  be  reproduced  to  properly 
estimate  storm  impacts,  since  field  surveys  may  be  made  after  some  recovery 
has  already  taken  place. 

221.  SBEACH  and  similar  beach  erosion  models  are  sensitive  to  changes 
in  water  level.  Considerable  effort  is  required  to  make  progress  in  predict- 
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ing  water  levels  at  the  coast  In  an  economical  and  complete  fashion.  Pheno¬ 
mena  to  be  studied  and  Incorporated  Into  a  coastal  engineering  methodology 
compatible  with  profile  change  modeling  Include  wind  setup,  infragravity  wave 
motion,  and  contributions  from  differentials  In  atmospheric  pressure. 

222.  In  the  present  version  of  SBEACH,  the  same  transport  formula  is 
used  Independent  of  transport  direction.  This  description  Is  probably  accep¬ 
table  In  the  Initial  stage  of  beach  recovery  when  waves  break  offshore  and  a 
well-established  surf  zone  exists.  However,  as  material  Is  transported 
onshore  and  waves  break  closer  to  shore ,  the  width  of  the  surf  zone  is  reduced 
and  che  character  of  wave  breaking  might  change  abruptly  from  plunging  to 
surging.  This  development  during  the  recovery  phase  is  not  described  satis¬ 
factorily  by  the  model.  Also,  onshore  transport  in  deeper  water  outside  the 
break  point,  which  takes  place  over  long  time  periods,  lacks  a  good  descrip¬ 
tion  in  the  model.  Finally,  regarding  the  transport  rate  formulation,  only 
two  types  of  cross -shore  distributions  are  presently  allowed,  purely  onshore 
or  purely  offshore  transport.  During  times  of  changing  wave  conditions  and 
mixed  waves,  more  complex  transport  rate  distributions  occur,  particularly  if 
the  profile  is  nearly  in  equilibrium  with  the  waves. 

223.  It  is  possible  to  specify  different  grain  sizes  (equilibrium 
energy  dissipation)  according  to  location  along  the  profile  in  the  model,  but 
since  no  tracking  of  individual  grains  is  performed,  the  description  is  rather 
crude.  Thus,  it  is  only  realistic  In  the  present  model  to  simulate  beach 
conditions  involving  a  few  distinct  grain  sizes.  By  keeping  track  of  the 
grain  movement  and  assigning  a  probability  distribution  to  the  grains  in  each 
calculation  cell,  it  would  be  possible  to  improve  this  aspect  of  the  model. 
However,  such  a  method  would  require  considerably  more  computer  memory. 
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APPENDIX  A:  NOTATION 


A  Shape  parameter  for  equilibrium  beach  profile,  ml's 
C  Wave  phase  speed,  a/sec 

Cg  Wave  group  speed,  m/sec 

Dj0  Median  grain  diameter,  mo 

0  Wave  energy  dissipation  per  unit  water  volume,  Nm/m3/sec 

E  Wave  energy  density,  Nm/m3 

F  Wave  energy  flux,  Nm/m/sec 

g  Acceleration  due  to  gravity,  m/sec3 
h  Water  depth,  a 

H  Wave  height,  m 

H  Mean  wave  height,  m 

K  Sand  transport  rate  coefficient,  a‘/N 

L  Wavelength,  m 

M  Dimensionless  empirical  coefficient 

n  Ratio  of  wave  group  speed  and  wave  phase  speed 

q  Net  cross -shore  sand  transport  rate,  m*/B/»ec 
S^  Radiation  stress  component  directed  onshore,  N/m 
t  Time,  sec 

T  Wave  period,  sec 

w  Sand  fall  speed,  m/sec 

x  Cross -shore  coordinate,  m 

y  Longshore  coordinate,  m 

ZR  Height  of  active  subaerial  profile,  m 
&  Local  beach  slope  seaward  of  break  point 
T  Ratio  between  wave  height  and  water  depth  at  breaking 
T  Empirical  stable  wave  height  coefficient 

A  Change  in  quantity 

e  Slope-related  sand  transport  rate  coefficient,  m3/sec 

•q  Wave  setup  or  setdown,  m 

it  Wave  decay  coefficient 

A  Spatial  decay  coefficient,  m"3 

6  Wave  angle  with  respect  to  bottom  contours 

Al 


> 


p  Density  of  water,  kg/m3 

4>  Local  profile  slope 

4>ly  Angle  of  initial  yield 

4>r,  Residual  angle  after  shearing 

C  Surf  similarity  parameter 

Subscripts  and  Superscripts: 

b  Breaking  condition 

eq  Equilibrium  condition 

es  End  of  surf  zone 

i  Variable  at  grid  point  or  cell  i 
k  Variable  at  time  step  k 

N  Number  of  cells  where  avalanching  occurs 

o  Deepwater  condition 

p  Plunge  point 

r  Runup  limit 

s  Stable  condition  for  depth- limited  wave 

1  Specific  value  of  variable  in  zone  I 

2  Specific  value  of  variable  in  zone  II 


A2 


